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Chapter 1 

Classical and Quantum Hall 
Effect 



The research in condensed matter physics is very often closely related to ad- 
vances in material science. The improvement of the experimental techniques 
to fabricate novel devices and the achievement of very low temperatures have 
opened new extraordinarily opportunities to investigate the physics of the 
electron systems. One of the most impressive result in the last twenty years 
is the discovery of high nonlinearities in the conductance of a two-dimensional 
electron liquid in the presence of a high magnetic field 0, Q, 0| • The name of 
Quantum Hall Effect (QHE) was assigned to this phenomenon because it re- 
calls, as we will briefly discuss later, the classical Hall effect in metals. This 
discovery has started a very productive research field where new concepts 
as Landau Levels, Incompressible Electron Liquid, Composite Fermions and 
Edge States were introduced and proved to be very useful 0, H, IE IE H 1^ ■ 
Also new experimental ideas were developed and nowdays the Quantum Hall 
Effect is used as a standard to define some physical constants. 



The experimental results for the transversal and longitudinal pxy re- 
sistance are summarized in Fig. 11.11 (data from [HI). The transversal resis- 
tance of this system shows an exact quantization in its dependence on the 
magnetic field, p^y = h/e^v with v an integer or rational. The longitudinal 
resistance p^x-, on the other hand, shows a very non trivial behavior, vanish- 
ing or having a delta-like behavior when the transversal conductance resides 
on a plateau or passes through two distinct plateaus, respectively, (see for 
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MAGNETIC FIELD (Tesia) 



Figure 1.1: The longitudinal and transversal resistance in Hall bar as func- 
tions of the magnetic field. In the transversal resistance the plateaus of 
quantized values are evident. From [l3 |. 



reference Fig. II. Ij) . The QHE can be expressed hence by the relations^ 



. eV / 1 \ (1.1) 



a 



h \ -1 



The first equation is the Ohm's law in matrix form and relates the current 
density to the electric field in the material. In the plateau regions the con- 
ductivity matrix becomes purely off-diagonal and is expressed in terms of a 
universal conductance quantum (e^/h) and of a number i>. From the experi- 
ments we know that u, which takes the name of filling fraction, assumes only 
integer and fractional values. We separate the case when u is an integer, and 
call this phase Integer Quantum Hall Effect (IQHE), or a fractional number, 
namely the Fractional Quantum Hall Effect (FQHE). 

It is widely believed that the IQHE and FQHE have a very different phys- 
ical origin. IQHE is thought to come from the interplay between the Landau 

^We will use the boldface to define vectors in two dimensions and an arrow to define 
three dimensional vectors. In the following we will also use a relativistic notation where 
the greek letters will indicate the coordinates in 2 + 1 dimensions and the latin letter the 
spatial coordinates. 
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Levels quantization of the motion of a charged particle in a magnetic field, 
and the disorder. FQHE cannot be explained without taking into account 
the electron-electron interaction. This leads to a new state of matter in 
which one may have a new kind of quasi-particles earring a fractional charge 
j^. A hierarchy of fractional filling fractions was justified theoretically by 
assuming that in the FQHE the quasi-particles which consistute the excita- 
tion of the system can undergo a IQHE. The idea is that every quasi-particle 
carries a number of magnetic flux quanta which partially compensates for 
the external magnetic field. Hence the quasi-particles may give rise to a 
IQHE phase when the effective magnetic field, obtained as the difference of 
the external applied field and the sum of the magnetic fields attached to the 
quasi-particles, assumes certain values j^. 

The quasi-particles defined in these theoretical approaches belong to the 
bulk of the system. However it was showed that the current in the QHE 
is concentrated on the edges of the system jlH. This can be justified by 
a semiclassical approach. The finite size of the sample can be taken into 
account by an external potential which bends up near the physical edges. The 
electrons that are well inside the system bulk move on circular orbits whose 
radius is determined by the magnetic field. Hence they are fully localized and 
cannot contribute to the transport. However the orbit of the electrons which 
reside near one edge is not fully contained in the system. From a classical 
point of view these electrons are reflected by the barrier generated by the 
confining potential and they propagate all in the same direction because the 
magnetic field fixes the sign of the circulation. Then there is the possibility of 
the formation of 'skipping orbits' which can travel through the system. From 
a quantal point of view the electrons localized near the edge of the system 
reside in the 'edge states' whose properties we will discuss in the following. 

In this thesis we will be concerned with transport properties and edge 
dynamics in the FQHE regime. We will show that the charged excitations 
of a two-dimensional electron liquid are concentrated in the edges of the 
systems thereby forming a one-dimensional system. Hence we can study 
the transport properties of the whole system by studying the interaction 
beetween two one-dimensional systems. It is believed that the edges of a 
system in the FQHE regime are a nearly perfect realization of the (Chiral) 
Luttinger Liquid (xLL)[l3, E^. The Luttinger Liquid model describes one- 
dimensional electron liquid with linear energy dispersion. It is an exactly 
solvable model and it will be a fundamental tool in this thesis. We will 
describe its basic properties in the following chapter. 

Starting from some physical reasonable assumptions we will show that the 
localization of the density excitations in the region of the edges is independent 
of the presence of a developed QH phase. We will derive an equation of 



5 



motion for such excitations and then use the solutions of this equation to 
calculate the electrical transport properties. 

This research was stimulated by a collaboration with an experimental 
group. The main experimental aim was the realization of new electronic 
devices able to detect the fractionally charged quasi-particles. To this end, 
it was planned to study the transport through a constriction in the two- 
dimensional electron liquid. This kind of problems have been treated the- 
oretically in the past by considering the constriction as a quantum point 
contact in the sense that in a very narrow region the edges of the Hall liq- 
uid come almost in contact so enhan cing the probability of the tunneling 
of excitations between the edges 0, 013. In this thesis we will drop 
the assumption that the constriction has zero extension, and we will show 
that, by taking into account a finite size of the constriction, one obtains 
different results for the tunneling current. In particular, we may obtain a 
better qualitative agreement with the experimental results in certain ranges 
of temperature. 



1.1 The Classical Hall effect 

We will start our discussion of the QHE by recalling the Classical effect. The 
Classical Hall effect occurs when a current flows in a metal in the presence 
of a uniform magnetic field. The motion of a charged particle in a magnetic 
field can be separated in the motion parallel to the direction of the magnetic 
field and in the motion perpendicular to such direction. We assume that the 
magnetic field is parallel to the direction of the z axis 

B = -Bz (1.2) 

and we choose a right handed set of orthogonal axis {x,y,z). The equation 
of motion projected in the z direction is given by 

z = -eE, (1.3) 

hence the particle is uniformly accelerated in this direction. 

We neglect in the following such motion and consider the particle confined 
in a two-dimensional plane {x, y) . The equation of motion can be written as^ 

mr = — eE — ev x B. (1.4) 

^The electron charge is negative then the physical constant e is positive. 
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When we consider the equihbrium situation we have that the acceleration 
must be zero and we obtain 

J = <tE (1.5) 
where we have defined the current density 

J = — nev, (1.6) 

and the conductivity matrix 

. f -aH\ 
"=U^ ) 

The classical Hall resistance (Rh), defined by 

<rH = Rh' = ^, (1.8) 

where n is the electron number density, is then linearly dependent on the 
magnetic field. 

The measure of the Hall resistance as a function of the magnetic field gives 
information on the electron density in metals. Notice that the Hall resistance 
is measured in the transversal direction with respect to the direction of the 
current. The conductivity matrix however is not complete. Indeed the Eqs. 
(jl.Sp and ()1.7|) imply that the longitudinal conductivity is zero. Hence the 
current in that direction will fiow without dissipation. In fact the classical 
Hall effect and the usual metal dissipation must be considered together then 
giving the conductivity matrix for a metal in a magnetic field 

a=f^^ -^^^ (1.9) 

where an is the Drude conductivity. 

We can understand the classical Hall effect in the following way. When 
the magnetic field is zero the electron current fiows following the gradient of 
the electric potential. If one turns on the magnetic field the motion of the 
particle becomes circular as one can easily verify be solving the equations of 
motion. The important point to stress is that the circulation is determined 
by the sign of the magnetic field and hence all the particles, which run in the 
same direction, turn in the same way. This effect creates a mean transversal 
current and a charge accumulation at one edge of the metal. To maintain the 
electrical neutrality an equal amount of opposite charges must accumulate at 
the other edge. This creates a transversal electrical field which will equilibrate 
the transversal current and establishes a dynamical equilibrium. The measure 
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(1.7) 



of the transversal electric field constitutes a way to measure the magnetic 
field. 

We introduce now a relativistic notation which will become useful in the 
following. The electromagnetic field can be defined by using the vectorial 
and electric potentials A and V. In terms of these potentials the electric and 
magnetic fields are defined as 

E = -dtA-VV, 

^ ^ (1.10) 
B = \/ X A. 

We deal with only two spatial dimensions because, as we have seen, the 
problem can be separated and in the direction of the magnetic field the 
equation of motion is easily solved. On the other hand when we will address 
the quantum problem we will see that the particles are strongly confinated in 
the direction of the magnetic field and the problem can again be separated. 
We introduce the covariant "tri-vector" potential^ 

A^ = {V,A) (1.11) 

where the index /i will run over 0,1,2 with the convention that will coincide 
with the temporal part of the tri-vectors, 1 coincides with the x variable and 
2 with y. As an example the tri-derivative is 

d, = {d,t, -d,, -dy) = {d,t, -V). (1.12) 

We need to introduce also the metric 



g^.u=\ -1 , (1.13) 




and the totally antisymmetric Levi-Civita tensor e^j/. With these definitions, 
we may rewrite the magnetic and electric field in the form 

E, = -doAi + d,Ao 
B = e.kd'A^ 

where, as usual in the relativistic notation, repeated indices are summed over. 
We can combine the equation ()1.5|) and the definition ()1.8|) in the compact 
form 

= ^He^^d-A"^ (1.15) 

•^Notice that even if the magnetic field points out the (x, y) plane the vector potential 
can be defined as a function only of the in-plane variables. 
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where = cp = —ecn. This equation is the phenonemological representan- 
tion of the Hall effect and it can be extended to the QHE because it comes 
directly from a Lorentz transformation which is valid if the system is transla- 
tionally invariant. Its quantum counterpart constitutes the starting point for 
the seminal papers of Wen , 0, 0, 0] for the derivation of the properties 
of the FQHE. The main task of a microscopic theory of the QHE must be to 
derive this equation starting from a microscopic electron Hamiltonian. 



1.2 A first glance to the Quantum Hall Effect 

The classical theory of the Hall effect was well understood and proved very 
useful in the experimental characterization of the properties of non-magnetic 
metals. Hence the experimental result of von Klitzing et. al. reported in 
Fig. 11.11 was totally unexpected. 

Before discussing the features of this experimental result, let us consider 
briefly the experimental setup. The devices used in this type of experiments 
are semiconductor heterostructures or hetero junctions where the electron 
gas resides at the interface between two different semiconductor species. The 
polarization of the different semiconductors generates a uniform electric field 
which confines the electrons in a small region localized around the interface 
as it is indicated in Fig. 11.21 Hence the electron gas can be considered as 



Al SiO Si 





conductance band 
3 nnri^— — ~ 


^^-"^^ valence band 


_l_ 









+ 



Figure 1.2: A simple scheme of the electron confinement around an interface 
between two different semiconductors. 



two-dimensional. This introduces a great simplification. Indeed it is well 
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known that the resistivity p, and the resistance R are related by 

R = pL^^-'i) (1.16) 

where d is the dimension of the system and in two dimensions these two 
quantities coincide. 

We choose the direction of the electric field as the z axis. In this direction 
a magnetic field will be also applied. To reduce the thermal effects the device 
is maintained at temperatures below 1 K (the newest experimental setups 
can reach about 10 mK). The electron gas is contacted with a multiprobe 
setup. Two of these contacts are used to inject a steady-state current and 
two other contacts are used to measure the potential drop in the system. A 
six-probe setup is shown in Fig. II. HI 




Figure 1.3: a) A schematic of the device used in the experiment on the Quan- 
tum Hall effect. The source (#1) and drain (#4) contacts are used to inject 
in the system a steady-state current /. The other probes (#2, #3, #5, #6) 
are used to measure the potential drop in the system to obtain the curves in 
Fig. O 

Once the current is fixed we can measure the Hall resistance by measuring 
the voltage drop, as for instance, on the probe 4 and 1 or the longitudinal 
resistance by measuring the potential of the contacts 4 and 3 of the device in 
Fig 11.31 The results of these measurements at very high magnetic field and 
very low temperature are shown in Fig. 11.11 In that figure the longitudinal 
and transversal resistances are reported. By using the Eq. ()1.8j) for the Hall 
resistance and the Drude part for the longitudinal resistance one expects 
that the longitudinal resistance is a constant and the transversal is linearly 
dependent on the magnetic field intensity. 

From the experimental data we see that this is not true (see Fig. 11.111 . 
The transversal resistance shows clear and large plateaus. In these plateaus 
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the resistance has a fixed value which is proportional to the quantum of 
conductance /h and the constant of proportionality is a integer or rational 
number. Outside the plateaus its relation with B becomes linear. More 
strange is the behavior of the longitudinal resistance. For small value of the 
magnetic field the Shubnikov-DeHaas oscillations are seen. However when 
the amplitude of these oscillations increases the longitudinal resistance shows 
an alternation of high maxima and fiat regions where its value is zero. The 
plot shown in Fig. 11.11 is a more recent version of the combined result of von 
Klitzing jl| on IQHE and Stormer, Gossard and Tsui [3] on FQHE. 

What is paradoxical in these results is that the electron-electron inter- 
actions and the disorder play a fundamental role to obtain it rather than 
destroy it. Indeed the results for the conductances in the previous section 
can be derived by using only the relativistic invariance jl^ , hence it does not 
matter if the system is classical or quantum and the electrons form a gas or 
a solid just maintaining the translational invariance. 

Let us notice another peculiar aspect of this effect. In the region where the 
longitudinal resistivity is nearly zero and the transversal one is quantized, 
the system behaves like a superconductor [p^x = 0) and like an insulator 
= 0). The reason for that is the presence of the finite off-diagonal 
terms, axy, which assures the existence of the inverse of the matrix a and the 
finite value of the pxy resistivity. 



Let us start the theoretical study of the QHE from the Integer Quantum 
Hall Effect, this being the first experimentally observed ^ and the more 
theoretically understood. In this effect the resistivity pxy is quantized in 
terms of the quantum h/e^ and the costant of proportionality is an integer 
(measured with a precision of a part over 10^). 

The very starting point is the quantum motion of a charged particle in 
a magnetic field. As it is well known, the motion of a charged particle in 
a magnetic field is obtained from the free Hamiltonian with the "minimal 
substitution" 



where A is the vector potential''. The presence of an electric potential does 
not affect such transformation and the Hamiltonian then reads 



*We choose to indicate the vector potential as a bidimensional vector. We assume that 
the component in the z direction is identically zero. 



1.3 The Integer Quantum Hall Effect 



p ^ p + eA 



(1.17) 
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We consider the simplest case of $ = 0. There are many ways to solve the 
Schrodinger equation for a single particle. We consider here a field-theoretical 
inclined approach that will provide some useful relations to understand the 
modern literature. The position and momentum operators follow the usual 
commutation relations 

[x,Px] = [y,Py] = ih, ^^^g^ 

others = 0. 



We define the new operators 



^ = ^ by + e^y) , 

V = -^(Px + eA^), (120) 

Y = y~r], 
which follow the commutation relations 

others = 0, 

where we have introduced the magnetic length i'^ = h/eB. To obtain these 
commutation relations we must remember that and Ay are in general 
function of both the operators x and y. On the other hand, the rotor of A 
is independent from x and y and hence the vector potential must be a linear 
function of these operators and the commutators turn out to be c-numbers. 
In terms of these new operators the Hamiltonian takes the simple form 

H='^u;l{e + V') (1.22) 

where we are faced with the problem of a one dimensional harmonic oscillator 
in the variables r] and ^ so that the energy spectrum is written immediately 
as 

E„ = (n + ^ ) fiu^ (1.23) 



where we have introduced the ciclotron frequency Uc = eB/m. The energy 
levels take the name of "Landau Levels" (LL) after the solution of Landau 
of the motion of a particle in a magnetic field [2o| . 

To obtain an explicit form for the eigenfunctions we must specify a gauge. 
A possible gauge is the Landau gauge 

A = {By,0) (1.24) 
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and we factorize the eigenfunction as 



il,{x,y) = -^Xiy) (1-25) 



where x is a solution of the Schrodinger equation 



2m dy'^ 2 



x{y)^Enx{y). (1-26) 



This is the Schrodinger equation of a harmonic oscillator with the center at 
kxi'^ hence the solution is readily written in terms of a product of a gaus- 
sian and the Hermite polynomials. Notice that when evaluated over these 
functions the mean value of Y depends only on 

{Y) = kj\ (1.27) 

and because Y does not depend on time we have that in this representation 
Y = kjH. 

Now let us turn on an electric field of amplitude Ey in the y direction. 
The electric field adds a term which is linear in y hence when we substitute 
y = r) + Y,we have a term linear in rj that changes the center of the harmonic 
oscillator and does not affect the energy, and a constant term in Y which 
enters the energy 

n+-]nu;,+ ^— I + eEy{Y). (1.28) 



From this relation it is easy to calculate the velocities defined as 



IdEn 



h 9k 

We have 



(1.29) 



■2 



_i_ 

"""" n^^' (1.30) 

Vy^O 



and we can calculate the conductivity 



J, = 
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where n is the electron density (A^ is the total number of electron) and we 
have defined the filling factor 

n N hn , . 

^ = — = ^ = 1-32 

ub Nb eB ^ ^ 

From the Eq. ()1.31|) one can directly compute the conductivity matrix as 

In the definition of the filling factor ()1.32j) we have also defined the quantity 
ub = eB /h = (27r£^)~^. It is possible to relate this density to the degeneracy 
of the Landau levels. In fact let us consider the case of an infinite system 
which is constituted by the replication of a finite size, let us say L^x Ly = S, 
system. The quantum number must be quantized in this finite size system 
as = 2Txj /Lx and, on the other hand, the center of the electron orbit must 
be contained in such a system hence we have the upper limit kx(? < Ly. We 
see then that the discrete number j, which defines the quantization of kx-, 
must be limited by 

e 

Notice that this degeneracy is the same for every Landau level. From this 
relation we see that the quantity Nb can be interpreted as the number of 
magnetic fiux quanta contained in the system (we have used that $o = h/e 
is the magnetic fiux quantum). Hence the filling factor can be interpreted as 
the mean number of fiux quanta carried by the single electron. On the other 
hand the filling factor gives information about how many Landau levels are 
filled, indeed every Landau level can contain, as maximum, j'^ electrons^ and 
directly from the definition we have v = N/jm- This simple model seems to 
reproduce most of the physics of the IQHE (the form of the conductivity 
matrix, its dependence of the universal quantity e^/h ...), however up to now 
the filling factor u is a positive real number and there is no way to introduce 
a quantization of this quantity. To introduce a quantization of the filling 
factor we need to consider an open system, i.e. a system where the number 
of electrons is not fixed but may vary due to the presence of one or more 
reservoirs. 

Before discussing this situation let us consider the case of the electron gas 
confined in a potential V{y) which is almost fiat at the center and bends up at 

^Notice that we consider the electrons fuhy polarized, hence there is not the factor 2 
due to the spin. 
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the edge of the system^. We also assume that this potential is smooth on the 
scale of the magnetic length. If we assume that the potential is translationally 
invariant in the x direction we can again separate the variables in Schrodinger 
equation and look for a solution of the form 



where fk{y) will not be the solution of the harmonic oscillator. Because of the 
hypothesis that the potential varies smoothly with respect to the magnetic 
length scale the function fk^ is centered around kxi"^ = Y^^ and the energy 
E will be again given by the kinetic energy plus a potential term given by 
V{Yj:^). Hence the group velocities are given by 



and again the Vy component will be zero. The energy will depend on the 
momentum kx only in the region of high variation of the potential V hence 
the particle will have a non-zero x velocity only if they are close to the edges 
of the system. This then defines two regions, the left and the right edge 
where the electrons move with opposite velocities (recall that the energy will 
decrease when entering the device, from the left, and increase when leaving 
the device to the right). The width of these regions depends on the variation 
of the potential: if the potential V varies very sharply these regions are very 
small, while if V varies smoothly they are very large. 

An example of this behavior is the simple case when the potential V{y) is 
infinite outside the physical edges of the device and constant inside (it may 
be also zero). Hence its gradient is concentrated only in two points. We are 
then faced with the problem of a particle in a high magnetic field confined in 
a finite region. In the Landau gauge this problem is mapped into the problem 
of a particle moving in a ID well and subjected to the harmonic potential. 
This problem is exactly solvable in terms of special functions (the Kummer 
functions M and U ) |2l| and the condition of the confinement gives a relation 
between the momentum and the energy En{kx)- On the other hand when 
the magnetic field is very large we can solve the problem in an approximate 
way. Indeed in this hypothesis the magnetic length i will be very small and 
when the electron center kxi'^ is far away from the boundaries the energy is 
simply given by the harmonic oscillator energy and the ground state energy is 
the corresponding lowest energy level of the harmonic oscillator. However if 
the electron center kxi"^ is exactly at one boundary the electron wave function 

^In fact the edges are created by the presence of this confining potential. 




(1.35) 
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must have a node in that point and be zero outside the system. This will 
imply that the ground state for this case is the first excited state of the 
harmonic oscillator. We then see that the energy increases continuously 
when the electron center moves from the bulk region, far from the edges, 
towards the boundary regions, close to the edges. A similar behavior is 
verified when the confining potential is smoothly varying. What happens 
when many Landau levels are filled? We can apply the above idea to every 
Landau level separately (we suppose that the electrons are not interacting) 
and the result is simply that there will be two edge states for each Landau 
level. 

Let us now consider the more realistic case of an open system, where two 
reservoirs with different chemical potentials can inject electrons in the gas. 
We choose to place the reservoirs along the x direction. What is now the 
total current flowing in the system? Because the currents at the edges flow 
in different directions we have simply 



j dii ev^g{ii) (1.37) 



where 

ldj{E) _ I dY 

is the one-dimensional density of states and we have used the fact that j is 
related to by = 27ij/L. We use the one-dimensional density of states 
because the electrons that are confinated at the edge are moving in one 
dimension only. By using the expression for the velocity one arrives directly 
to 

^ = ^(/UL-/iii). (1.39) 

This current is the same for every Landau level hence the total current will 
be 

eN e^N 
It = NI = —{fiL - /ii?) = -f^{VL - Vr) (1.40) 

where Vi is the potential of the i-th. edge (/ij = eVi). The integers N count the 
number of Landau levels that are beyond the Fermi energy (we are assuming 
that the chemical potentials hl and are close enough to consider only 
the linear response regime). Now we can calculate the conductivity matrix. 
When the probes belong to the same edge the current is the same (there is 
not backscattering) and one obtains a zero longitudinal resistance. When 
the probes belong to different edges the current is given by p.4()|l . We then 
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obtain 




(1.41) 



It is now possible to understand the nature of the plateaus. The reservoirs 
will define an energy reference (the Fermi energy) and moving the magnetic 
field will vary the frequency ujc and the Landau level energy. Only the Landau 
levels that are below the Fermi energy can participate to the current and this 
defines the number N . When a new Landau level sinks below (or become 
greater than) the Fermi energy the number N will increase (decrease) by 1 
and the conductance can vary only by integers. On the other hand there is 
a finite energy gap, given by tiWc between two different Landau Levels. This 
finite gap implies the finiteness extension of the Hall conductance plateaus. 

What remains unexplained by this model is why when a new Landau level 
sinks below the Fermi energy the longitudinal conductance has a maximum 
and the Hall conductance shows a linear behavior. It is mostly believed that 
this behavior can be explained by inserting in the above model the effect of 
the disorder. The idea is that when a new Landau level crosses the Fermi 
level, a path which connects two edges may appear and this allows for the 
presence of backscattering and the longitudinal conductance is not zero. The 
connection between the two edges is created by the presence of the impurities 
that create localized islands and the electrons can, tunneling from one island 
to the other, be backscattered to the other edge. Numerical evidences of 
this percolation scheme have been provided but there is not an accepted 
model which reproduces the experimental data and allows us to understand 
completely this phenomenon. 



1.4 The Fractional Quantum Hall Effect 

The theory we have presented in the preceding section can explain the IQHE. 
However in the 1982 Tsui et al. showed that in cleaner sample and at higher 
magnetic field than that used by von Klitzing a new series of plateaus in 
the Hall resistance appears 0](see Fig. II. ip . These new plateaus appear 
when all the electrons are in the lowest Landau level, hence there is not the 
possibility that new edge states can be created when varying the magnetic 
field. Another obscure point is that the transversal conductance is again 
given by an expression similar to Eq. 1)1.411) but with the integer number 
substituted by a rational number. This rational number is identified with 
the filling fraction v. This effect is now widely known as Fractional Quantum 
Hall Effect. 
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There is not a complete accepted theory for this effect even though many 
theoretical ideas have been developed to understand its physics. A review 
of all these ideas is outside the scope of this thesis. However the interested 
reader can refer to some beautiful reviews [13, HI, 2^ > ^ • 

The first observation we can make in dealing with this effect is that the 
kinetic energy, since all the electrons reside in the LLL, is a constant and can 
be neglected. Hence the electron-electron interaction, which in the IQHE is 
treated as a perturbation, plays a fundamental role. If the interaction cannot 
be neglected the correlation between the particles must be taken into account 
in the search for the ground state. When we consider the problem of many 
electrons usually we first write the single particle wave function and then the 
many-particle wave function is given by a Slater determinant. However when 
the particles are strongly interacting this approach fails and we must write 
a wave-function which takes into account the correlation. On this track, 
Laughlin moved in the 1983 [3]. He wrote a trial wave function of the form 



|2 



V^ = n(^.-%rexp , (1.42) 

where m is an odd integer^ and Zj = xj + iyj is a complex variable which 
describes the position of the j — th particle. He proved by using numerical 
calculation that this wave-function minimizes the interaction energy. This 
minimization can be easily understood by observing that in the Laughlin 
wave-function the multiplicity of the zeros is greater than that of the fermion 
(this case is recovered by setting m = 1). Laughlin was also able to show 
that the state described by this wave-function is incompressible, hence there 
is a gap in the energy spectrum. The presence of the gap is necessary to 
understand the finite extension of the plateaus as in the case of the IQHE. 

Laughlin pointed out also that the quasi-particle excitations of this ground 
state have fractional charges given by e/m. The direct observation of such 
fractional charges has captured the major experimental efforts and was achie- 
ved in the 1997 by several experimental groups with different techniques 

The Laughlin's proposal can explain the presence of the plateaus at filling 
fraction given by a rational number which is the inverse of an odd number. 
However a more recent theory extended the model capturing also a series of 
possible values of the filling fractions. The starting point of these theories 
was the pioneering work of Jain which pointed out that the FQHE can be 



^The request that m is an odd integer is due to the anti-symmetry related to the Pauli 
exclusion principle. 
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explained as the IQHE phase of the ground-state excitations. These quasi- 
particles, called Composite Fermions, carry a certain number of magnetic 
flux quanta per particle and see a residual magnetic field. Under certain 
conditions these fermions can develop a IQHE phase. With this idea Jain 
was able to predict a series of filling fraction given by 

V = — - — (1.43) 
pm ± 1 

where p is an odd integer and m is an integer. 

An earlier approach to the generalization of the allowed values of u is 



the hierarchical approach developed by Haldane j26[ and Halperin j27|. The 
general idea of this approach goes as follows. We consider a given FQHE 
state in its ground state. When we change the magnetic field we add quasi- 
particles to the system. These quasi-particles move in the magnetic field 
determined by the difference of the external magnetic field and the magnetic 
field carried by the particles which form the ground state. When we increase 
again the magnetic field the quasi-particles develop a new Quantum Hall 
phase. The new filling fraction is now given by a complex fraction composed 
by the integer numbers which identify the starting state and the new integer 
which identify the new developed phase. A complete review of these ideas 
can be found in Wen A major aspect of this theory is that a series of 
excitations starts to develop and a given filling fraction u = m/p might be 
the result of many of these steps. This in turn implies that many modes can 
propagate in the system. This can be the basis to understand some recent 
experimental results [1^ H, 0] 

A major shortcoming of this theory is that its starting point is the phe- 
nomenological relation ()1.15|) (see Ref. and references therein) where the 
presence of a developed Quantum Hall phase has been assumed implicitely. 
Up to now, at the best of our knowledge, no one was able to derive this equa- 
tion from the fundamental properties of the electrons in the magnetic field 
confined in a finite size device. Even in the Composite Fermions approach, 
one of the most widely used theories to study the FQHE, it is not clear how 
to derive, starting from the correlated wave function, the presence of a gap 
which can justify the quantization of the Hall conductance. 



1.5 Edge dynamics 

In a series of beautiful papers Wen shows that the edge dynamics can be 
described by using a so called Chiral Luttinger Liquid (xLL) [l^, 17, 14 1. 
We want to discuss briefly this theory without entering in the details. In the 
following chapters we will compare this theory with our model. 
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The Wen's starting point is the phenomenological equation (jl.lSp . This 
equation is the requirement that a QHE is fully developed and we are in 
one of the plateau of the Hall conductance. It relates the response of the 
system (the current J'^) to the external perturbation (the potential A'^) in a 
non-trivial way. We can define a Chern-Simons field a'^ by 

J'' = ^ne^uxd'a^ (1-44) 
and obtain the Eq. fjl.l5|l as the equation of motion of the density lagrangian 

C = —e^'^^A.dxAr,. (1.45) 

where A^ = a^j^ — A^ and is the dynamical field. 

As one can easily verify this lagrangian is not invariant with respect to 
a gauge transformation in a space with boundaries. To make the action 
gauge invariant Wen assumes that the lagrangian is written for the bulk 
operators and does not take into account the effects of the edges. Taking into 
account this contribution adds a term whose form can be derived from the 
gauge invariance and can be expressed in terms of current-current correlation 
functions. Wen was also able to show by assuming the locality of the theory 
that the excitations must have a gapless and linear spectrum. In particular, 
if one defines the current in the momentum space 

= I da^e'^'^J^ia), (1.46) 

where a = (0, 1) indicates the temporal (0) or spatial (1) component of the 
current vector and cr parameterizes the boundary, then one gets 

[H,J^]=ckJ^, 

["^k ^ "^k'] ~ ["^k ' "^k'] ~ 0' 

where by definition = 1/2(J° ± /c). In these equations we have as- 
sumed that c is the velocity of the boundary excitations and that they are 
chiral waves. From these equations it is also easy to show that the density 
fluctuations localized at the edge follow the commutation rules 

[4°,4] = — W,o. (1.48) 

Wen now shows that one can "bosonize" the theory. The bosonization tech- 
nique is a powerful idea to solve the problem of interacting one- dimensional 
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electrons model or Tomonaga-Luttinger model. The main observation is that 
in the low-energy regime, the densities of a one dimensional electron system 
are boson operators. Hence instead of solving the problem for many inter- 
acting electrons with fermion properties one can solve the problem of the 
density modes which have bosonic properties [l^ . 
We define a boson field with density lagrangian 

C = \[{d,cl^f-{d^ct>f] (1.49) 

and we require that this field is chiral, i.e. it satisfies the constraint 

(9o - d^)ct> = 0. (1.50) 
The equation of motion for this field is 

{d^-d,){d^ + d,)<p = Q (1.51) 

which can be solved by assuming that 
= ci>o + ^o+P^{t + a)+p^{t~a) + Y, \ Ke-'"(*-'^) + a„e-^"(*+'^)) (1.52) 

where 00, 00, P</)5 and define the zero wavelength mode. The conjugate 
momentum 

^ = 7^ = ^00 = P<A+P<^+y] («ne-'"(*-'^) + «„e-^"(*+'^)) = ct>L+ct>R (1.53) 

must satisfy the commutation rule [0, tt] = 1 hence we can derive the com- 
mutation rules of the operators in the expansion for 

[0o,P</.] = [0o,P0] = ^, (1-54) 
others = 0. 

It is possible to use the chirality condition to eliminate the field 0, pcf, and 
Cin- The hamiltonian is diagonalized 



H =^pl + Y,'^i^n (1.55) 

n>0 

where we have neglected the constant terms. 
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By defining 

r = ^e^^dp(t>L, (1.56) 

where 0^ denotes the modes which satisfy the chirality condition, we can 
verify, starting from the boson algebra, that these operators satisfy the alge- 
bra of the edge excitation. Hence we have obtained a boson representation 
for the density fluctuations operators. To fully define all the algebra we need 
also the operators for the quasi-particles. Indeed the density J" describes 
the edge density fluctuation with respect to an equilibrium density po- We 
look for an operator that can change such equilibrium density. To obtain a 
form for this operator we define the equilibrium charge 



g = ey dar{a) (1.57) 
and the quasi-particle operator 

[^,g]=e*^ (1.58) 

where e* can be different fron the electron charge. In fact in a physical 
system the quantum charge is the electron charge. However, here, we are 
postulating that the quasi-particle charge may not be the electron charge. 
We have seen in the preceding section that the "charge fractionalization" is 
useful to understand the phenomenology of the FQHE. 

In the chiral boson theory the charged operators have a form 

m =: e'^'^^ : (1.59) 

where the symbol : . . . : denotes the normal ordering and 7 is a constant. 
This constant can be determined by the commutation relation ()1.58jl by using 
the boson commutators p.54|l obtaining 

e* 1 

7 = ^- (1-60) 



Require that the operator \I'''" creates a unitary charge corresponds to choose 
= 1. The other requirememt that one can make is that the operator 
\I' satisfies the usual fermion (anti-)commutation relations. By using the 
Haussdorf lemma 

e^e^ = et^'^le^e^ (1.61) 

we obtain 7^ = l/u = odd integer. Notice that this last relation says to 
us that if the condition u = 1/odd integer is not fulfilled hence the fermion 
creation operator must be composed by many different branches. 
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One can now calculate the fermion correlation function G defined as 



G=(vl>t(t,a)vl>(0,0)). (1.62) 

By using again the Haussdorf lemma we obtain, in the limit of large system 
size, a/L 

1.6 Transport through a constriction 

It is possible to use the formulation for the edge dynamics we have depicted 
in the preceding sections to understand the effect of tunneling between the 



edges [1J|. Such a phenomenon is present when two edges are close enough 
to have a non zero superposition of the quasi-particle wavefunction. The 
inclusion of tunneling in the edge model is done "by hand" in the sense that 
there is not a completely accepted model to derive the tunneling amplitude 
r from the electron or quasi-particle properties. In the hamiltonian for the 
edges we then add a term proportional to the product of two distinct quasi- 
particle operators 

Ht = T^<\^<r + T*¥^^!l. (1.64) 

In the framework of the linear response theory one can calculate the tunneling 
current and the tunneling conductance which changes the Hall conductance 
in a non-universal way. The result of this approach is a non linear behavior of 
the current as a function of the tunneling voltage, the bias difference between 
the edges, i.e. the Hall voltage, and of the temperature In particular 

Wen obtained that 

I oc V^"-^ (1.65) 

where a = 2/ v. A similar relation holds for the temperature dependence of 
the tunneling current 

IocT'^-\ (1.66) 

These relations stem directly from the power-law behavior of the correlation 
function G of the quasi-particle. 

The non-linear relation between the current and tunneling voltage are 
expected also for the Luttinger liquid model However the intrinsic 

difficulties in obtaining a very clean one- dimensional electron liquid reduce 
the possibility of the experimental observation of these results. In this sense 
the Wen's papers open the possibility of a verification of these result on the 
edge of a quantum Hall liquid which, as we have briefiy discussed, forms a 
very clean one-dimensional Luttinger liquid. 
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Similar results were extended later by Kane and Fisher, by using a Renor- 
malization Group (RG ) approach, to the effect of the impurity which medi- 
ates the tunneling [l^ . In particular, they showed that any impurity between 
two edges under the RG analysis implies a strong coupling i.e. the edges will 
eventually close, the tunneling current diverges and the Hall bar separates in 
two distinct regions (this phenomenon is known as "overlap catastrophe"). 

In this thesis we will generalize this approach in many ways. The first 
point we want to make is the generalization of the approach to the edge dy- 
namics eliminating the limitation of incompressibility of the Hall liquid in 
the bulk. This is done by a hydrodynamical approach to the density fluctu- 
ation in the QHE and by the projection in the Lowest Landau Level. The 
projection on the LLL simplifies the Hamiltonian for the system quenching 
the kinetic energy but, on the other hand, introduces a non trivial quantum 
commutator between the density fluctuations and restores the full hardness 
of the problem. The hydrodynamical approach eliminates the fluctuations 
present at the scale of the magnetic length. This allows us to fully describe 
the system by using only the density fluctuations neglecting the real electrons 
and allows us to consider only the low energy excitations of the system. Sim- 
ilar approaches were used to study the tunneling of electron from a Fermi 
liquid to the edge of a fractional Quantum Hall liquid 2^ 3^ • 



On the other hand there are not conceptual difficulties in inserting in 
our model the intra- and inter-edge interactions. In this way we arrive to a 
problem similar to a Luttinger liquid model where an edge maps to a chiral 
branch of the liquid. However our edges are still interacting while in the 
Luttinger liquid model the reduction to chiral waves is done by eliminating 
the interactions. 

The scheme we have depicted before and that we will discuss in the fol- 
lowing brings us to the Hamiltonian for the density fluctuations. We bosonize 
these fields and obtain an equation of motion for the amplitude of the fluctu- 
ations. We choose a situable form of the intra- and inter-edge interaction and 
solve the equation of motion in many different edges profiles. In particular 
we are interested to the case when a constriction, which brings the edges to 
stay close, is present. It is possible to solve this case by using a scattering 
approach, i.e by defining the transmission and reflection coefficients for the 
wave which impinges on the constriction. We will show that the constriction 
does not change the low frequency (long wavelength) behavior of the con- 
ductance thus recovering the linear relation between the Hall conductance 
and the filling factor^. To fully understand some new experimental results 



^Notice that in our model the filhng factor is not quantized at ah, hence we are unable 
to recover the quantized Hall conductance. 
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j28| we will introduce the tunneling between the edges. We restrict the tun- 
neling only to the region of the constriction, hence we use the quasi-particle 
operators that are given by the solution of the equation of motion when the 
constriction is present. We develop a perturbative approach starting from 
the equation of motion for our operators and use it to calculate the tunneling 
conductance. This quantity, which modifies the Hall conductance, is a non 
linear function of the tunneling voltage and of the temperature as shown by 
WenQ. However the constriction introduces an energy scale. The charac- 
teristic time to is given by the time that an edge excitation needs to travel 
through the constriction. This time introduces the energy scale h/to. We 
will show that the response functions of the system have different behaviors 
if the energy is above or below this energy threshold and this will affect the 
tunneling resistance. 

Recent experimental results, reported by Roddaro et al. [2^, show that 
the tunneling conductance for low bias voltage has a different behavior when 
compared with the Wen's results. In particular they report that in their 
experimental setup at relatively high temperature the conductance shows a 
large maximum for zero bias voltage and two deep minima when the voltage 
is increased. When the temperature is lowered two main aspects appear. 
The zero bias maximum initially start to increase as expected but beyond a 
temperature near about 400 mK the conductance starts to develop a large 
zero bias minimum. Beyond this temperature also a strong asymmetry starts 
to appear. From the discussion with the experimentalist we known that 
the appearance of the central deep in the conductance is not related to the 
presence of localized impurities. 

Our model seems to capture some of the aspects of these results. In par- 
ticular the presence of the deep minima seems to be related to the presence 
of the constriction and the comparison with the experiments can give infor- 
mation on to and e*. We predict also the presence of small oscillations on 
the tunneling resistance at large bias. The frequency of these oscillations is 
given by 

e*to e*d cosh 26^2 

where Ci is the velocity of the modes out of the constriction and 9i and 62 
are related to the intra- and inter- edge interactions inside and outside the 
constriction. 

The thesis is organized as follows. In this first chapter we have discussed 
the general phenomenology of the QHE and the Luttinger Liquid model. To 
do that we moved from the Classical Hall effect, given a look to the main 
experimental results of von Klitzing Q and T^ J. Stormer, and Gos^ard 

and discussed the theoretical result for the IQHE and the FQHE. We 
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have also briefly reviewed the Wen's theory for the edge dynamics and then 
discussed our main results. 

In the second chapter we will concentrate on the derivation of our model 
and of the equation of motions. We will discuss its properties in the Hilbert 
space and solve the equation of motion in some simple cases that will prove 
to be useful in the following. 

In the third chapter we will discuss the problem of defining and calcu- 
lating the transport properties in these devices. We give a definition of the 
conductance and show how we recover, in the simple cases studied in the 
second chapter, the ideal Quantum Hall conductance. 

In the fourth chapter we will abandon the limitations of a translational 
invariant system and discuss the presence of the tunneling between the edges. 
In this case we show that the tunneling gives a non-universal correction to the 
Quantum Hall conductance. We also compare our model with those present 
in the literature. 

Finally, in the last chapter, we will discuss the comparison with some ex- 
perimental results and the future perspectives of this research, and sumarize 
our main conclusions. 

The results presented in this thesis have been submitted for pubblication 
in the Physical Review B journal. A preprint of this paper is available [8ll |. 
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Chapter 2 

The model for the Quantum 
Hall bar 



We are interested in the study of a two-dimensional electron gas confined in 
a finite region in the presence of a high magnetic field. We will show that the 
description of the electron gas in terms of two interacting chiral Luttinger 
Liquid follows from some semi-classical assumptions we will discuss in the 
following. We want to point out also that we do not require that the system 
is in a Quantum Hall phase. 

The ground state of the electron gas is characterized by a density po(i") 
which takes into account the electron-electron interaction and the action of 
the external potential which confines the system. We consider the excitations 
of this system as density fluctuations Sp{r) near the "equilibrium" density 
po{r). The dynamics of these fluctuations is determined by the Hamiltonian 



and by the commutation relations for the fields 5p. In considering the Hamil- 
tonian ()2.1|1 we have projected all the operators in the lowest Landau level. 
Such a projection quenches the kinetic energy (it becomes simply a constant) 
but, as we will discuss later, introduces non-trivial commutation relations for 
the field 6p. 

In the definition of the Hamiltonian (see eq. ()2.H) ) we have inserted an 
interaction potential V{r — r') between two density fluctuations at different 
points. This electron-electron potential has a coulomb origin. In the follow- 
ing, however, we will assume a simple form for this potential. 




(2.1) 
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2.1 The Luttinger Liquid model 



In this section we will briefly review the solution of the Luttinger Liquid 
model because we will use a similar approach to the solution of our model 
and this analysis will be useful when we will compare our model with others 
present in the literature. We do not derive the Luttinger model starting 
from the physical properties of the interacting electrons in one dimension. 
The interested reader can refer to the seminal paper of Haldane ^2] or to 
the book of Mahan [ij. 

The Luttinger Liquid model describes the dynamics of a one dimensional 
electron liquid. The Fermi surface of such a system is composed by two 
points then with the same energy there are two species of excitations, 
the left and right movers deflned by the sign of their momentum. We make 
the association Left = L = +1 and Right = R = —1. In the following we 
will consider also a spinorial notation where the upper component will be 
always related to the left movers. 

When the electrons are not interacting the Hamiltonian operator is (we 
use h = c = 1) 

Ho = vfJ dx{pl + pl) (2.2) 

where vp is the Fermi velocity and {pi) is the density operator of the 
right (left) movers 0, ^. The density operators follow the commutation 
relations 

[pa{x),pp{x')] = -ialpdj{x - x'), (2.3) 

where a and (3 assume the values R or L. These relations can be derived by 
starting from the deflnition of the density operators in terms of the electron 
creation and annihilation operators and their anti-commutation rules. 
It is customary to deflne the operators 

dJa = -apc,. (2.4) 

These new operators will be useful in the following, when we will introduce 
the interaction, to diagonalize the total Hamiltonian. It easy to show that 
Pa and 9a are conjugate fleld 

[Paix), Of^ix')] = -i5a,fs6ix - x'). (2.5) 

^This Hamiltonian differs by a constant from the usual definition one can find in the 
references. 
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When we consider the interaction we add to the free Hamiltonian Hq the 
term 



Hi =Vi / dx{p% + pI) + V2 dx{pRpL + PlPr) 

\ ^ r (2.6) 



where Vi and V2 are the hmit of vanishing momenta of the Coulomb in- 
teraction between the right and left movers. Vi is the interaction between 
electrons that belong to the same species, while V2 is the inter-species inter- 
action. Notice that the interaction preserves the momentum. 

The first step towards the solution of this problem is the diagonalization 
of the total Hamiltonian H — Hq + Hi. We can accomphsh this task with 
the linear transformation 



On = Or + 9l, 
Oj = Or — Ol, 



(2.7) 



and obtain the new Hamiltonian as 



H=^-^ I dx 



where we have defined 
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^-^-^/ti + -) -(-) , 



(2.9) 



' vf + V1 + V2 
^~ VF + V1-V2 



The Hamiltonian is separated but the new field ^jv and 9j are conjugate 
fields, indeed their commutation relations are 

Wnix), ej{x')] = [djj{x), eN{x')] = -2t6{x - x'). (2.10) 

Notice also that these fields are not chiral modes. The decomposition of the 
problem in terms of non-interacting chiral modes is not yet complete. 

To obtain a set of independent fields we consider the new transformation 
on the fields On and Oj 

= qOn + Oj, 
(t>L = qOn - Oj. 
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The Hamiltonian in terms of these fields is now 

H = ^Jdx [{dM' + (2.12) 

and we have 

[d,Mx),M^')] = -42(7<^5(x-x'). (2.13) 
The equation of motion for the fields are 



(2.14) 



hence these modes are chiral. 

We finally can express the density operator in terms of these chiral modes 



(2.15) 



f \ , g + 1 
P^^^' " ""^^ 4^ ""^^^ 

and obtain the total density operator p = + as 

P = Id^R - (Pl) = ^dt{<PR + <Pl). (2.16) 

The total current is defined by using the continuity equation p = —d^I and 
then 

/ = -|5.(<^i? + <^l) = -ldti(j)R - (j)L). (2.17) 

When we consider the chiral Luttinger liquid (xLL) we fix one of the 
fields to zero and consider the dynamics of the other field given by the 
Hamiltonian ()2.12j) and the commutation relation ()2.13p . Notice that we 
have started with two interacting chiral modes and after a non-canonical 
transformation we have obtained two non-interacting chiral modes with non- 
trivial commutation relation. As it is seen from Eq. ()2.15p . the chiral modes 
(J) are a linear combination of the starting interacting chiral modes 6r and 



2.2 Lowest Landau level projection and the 
hydro dynamical approximation 

The Hamiltonian ()2.H1 alone cannot determine the dynamics of the fields 5p. 
To do that we need to give either the equation of motion or the commutation 
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rules for these fields. To calculate the commutation relations we project all 
our field on to the lowest Landau Level. The theoretical approach to this 
projection and its physical implications are discussed in a more detailed way 
in many reviews or articles (see as an example Ref. [1^). In the following 
we will use it in fixing the wave functions of the single electron. 

In order to emphasize the role of the two main physical approximations, 
i.e., the LLL projection and the hydrodynamical approximation, it is con- 
venient to start from the second quantized form of the density fluctuation 

5p{v) = Y,Ach^l{^)Vh{v) (2.18) 

k^h 

from which we have 

[5p(r),5p(r')] = ^ {c\ci5h,rn - c{crn5k,i^ 

k^h,m^l (2.19) 

The indices k and h in Eq. fj2.18|) label states in the LLL within the Landau 
gauge A = (i??/,0,0). The condition k ^ h excludes the equilibrium contri- 
bution to the density po(r)- In the spirit of the hydrodynamic approach, the 
latter is related to the ground state expectation value 

4^/ {clck)5k,i = n{k)5k,i (2.20) 



via 

CO 

iy + kl')'^ 



f dk nik) 



/2 



(2.21) 



-oo 



Here n{k) is the occupation number for the state with momentum k and in 
the homogeneous case, n{k) = hq, the evaluation of the gaussian integral 
gives 

no 
2ttP 

where ly is the filling factor. 

The commutator between two different density fluctuations can be derived 
from the well known result 



Po = Tdn ^ '^o = 27r/Vo = (2.22) 



[p(q), p(k)] = (e^*^^^/2 - e-'"*''/') p(k + q). (2.23) 

When we consider the long wavelength density fluctuation we expand to the 
leading order in ki and qi and transform to real space to obtain 

[p(r),p(r')] - tehi,diPo{r)d,6{r - r'). (2.24) 
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Notice that we have substituted the density operator p with its equihbrium 
expectation value Po(r)- This approximation is justified as long as we are 
interested to the linear dynamics of small fluctuations about the equilibrium 
state. 

From these commutation relations and from the Hamiltonian ()2.H) we can 
derive the equation of motion of the density fluctuation 



By this equation it is evident that the density fluctuation are localized in the 
region of variation of the equilibrium density po- The regions of maximum 
variation of po are located near the edges hence the density fluctuation are 
conflned near the edges of the Hall bar. Because this model for the edge 
dynamics agrees with that we flnd considering the large magnetic fleld limit 
of the hydrodynamical Euler equations 34] we call our approach "hydrody- 
namical" . 

Attached to a point of every edge we consider a local right handed system 
of coordinates where the variable y indicates the region of variation of the 
equilibrium density. To be more precise we consider the device reported in the 
Fig. 12.11 In this flgure we draw a simple schematization of a Quantum Hall 
bar connected to two reservoirs (the drain (D) and the source (S)) bar where 
the physical edges are considered as parallels and translationally invariants. 
The coordinate y is then orthogonal to the edge while the x coordinate run 
along the edge. 

The observation about the localization of the density fluctuation allows 
us to greatly simplify the problem. Indeed rather than solve the complete 
two-dimensional problem determined by the equation of motion ()2.25p . which 
implies the exact knowledge of the equilibrium density po SI, we integrate 
over the y direction, starting from the edge to well inside the bulk when we 
consider the fluctuation localized on the left edge and viceversa for the right 
edge, and obtain a one dimensional density fluctuation. In this way we deflne 
the operators 



where the index a identifles the edge to which the density fluctuation be- 
longs. These new operators must satisfy the commutator law ()2.24|) when we 
integrate with respect to the y and y' variables, following the prescription in 
()2.2fi|l . It is possible to show that the result of this operation is 




(2.25) 




(2.26) 




(2.27) 
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Figure 2.1: a) The scheme of a simple Quantum Hall bar contacted with two 
reservoirs. If we neglect the effect of the disorder, the edges runs parallels 
and translationally invariant along the x direction. The y axis is chosen to be 
orthogonal to the edges, b) A magnification of the edge region. We integrate 
with respect the y direction to obtain a field varying only in the x direction. 



In the appendix El we give an alternative derivation of this equation starting 
directly from the operators 6pa in the case of sharp and smooth edges. We 
have also allowed for the possibility that the filling fraction u is dependent 
on the position x along the edge^. In the following, for simplicity, we will 
drop this degree of freedom and consider a uniform filling fraction. 

The Hamiltonian for the operators 6pa is derived from the Hamiltonian 
()2.H) by an integration. To do that we suppose that the interaction potential 
is slowly varying on the dimension of the edge i.e. we approximate 

V{r-r')c:.V^,p{x,x') (2.28) 

with 

V^,pix,x') = Vix,y^;x\yf,) (2.29) 

where we have allowed for the possibility that the two density fluctuations 
belong to different edges. With this approximation we get (repeated indexes 
are summed over) 

H=-J dxdx' 5pa{x)Vafi{x,x')5pp{x'). (2.30) 

^This can account, as an example, for the situation where two HaU hquids with different 
filhng factors, are connected. 
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This Hamiltonian together with the commutation relations ()2.27|) determines 
the equation of motion for the 5pa operator 

oo 

dtSpa{x) = — ^ / dx'dxVp^y{x,x')6py{x'). (2.31) 

Zn J 

—oo 

This is our first major result. The solution of this equation describes the 
excitations of our system and will be the fundamental tool to obtain the 
response functions. 

The commutation relations ()2.27|1 when u is constant are the reproduc- 
tion of the Kac-Moody current algebra ()2.13|1 for the density field in the 
xLL model. Notice however that our derivation of the xLL model has noth- 
ing to do with the QHE. Indeed it depends only on the high magnetic field 
limit (projection on the LLL) and on the coarse-graining (hydrodynamical 
approximation). Our derivation is therefore valid for every real value of the 
filling fraction whereas the QHE occurs only for certain value of u. We want 
to point out also that in our model every edge exhibits a xLL behavior even if 
the interactions are turned off, and that it is not possible to map the problem 
of two interacting QH edges to one Luttinger Liquid model. 

2.3 The equation of motion 

Let us now discuss the solutions of the equation of motion ()2.3H1 . To do that 
it is convenient to define the field (pa such that 

6paix,t) = d^(f)aix,t). (2.32) 

These fields satisfy the commutation relations 



(x), (j)p{x')] = i£-a^^^sign(x - x'). (2.33) 



In the following, we will consider also the function ifa{x) which are the solu- 
tions of the equation 



oo 

V 



iua^^pdxipisix) = — J dx' dxVa,i3{x,x')d^'ipp{x') (2.34) 

— oo 

and we will show that the field 0^ can be expressed in terms of the functions 
(fia and the operators b which turn out to have a Bose statistics. 

To complete this task, we need to study the properties of the solutions of 
the equation ()2.34j) . The first point we want to make is that these solutions 
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form a complete and orthonormal basis of the Hilbert space. More precisely 
we will show that: 

• all the frequencies ujn are real and we label these states such that uon > 
if and only if n > 0; 

• the eigenf unctions (pn,a form a complete basis in the Hilbert space with 
the completeness relation 

sign{uJn)<^nA^)VnA^')^^'p = " (2-35) 

n 

and the orthonormality condition 

Y^a,l3 / dx iflAA^x'^mAA = ^'^&^{'^n)^n,m\ (2.36) 

• the eigenfunction f^nAA is doubly degenerate. 

The proof of these properties is detailed in the appendix |Bl 

Having defined a basis for the Hilbert space we can develop the field 0^ 
and 5 Pa in this basis. Indeed for every (spinorial) function fa{x) we have 



fa{x) = Y dx' 6pA{^ - ^')fli.^') 



(2.37) 



where we have defined 

/(n) = -isign(u;„) J dx' v3n,a(a;'X,/3<9x'//3(a;')- (2.38) 
This result implies that we can define the operator 0(n) as 

5p„(x) = ^0(n)v9n,a(a;), (2.39) 

n 

where 0(n) turns out to satisfy the commutation relation 

[<i){n),<i)\n')] = ^sign(w„)5„,„., (2.40) 

ZTT 

and this suggests to define the operators b such that 



= {bnein) + blji-n)) . (2.41) 
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It is easy to show that we have [&]i,&n] = 1 thus these operators follow the 
Bose statistics. We have then proved that we can expand the 5pa operator 
as 

^Po.{x) = J^Yl (bnd.^nA^) + bld^^ijx)) (2.42) 
and we obtain for the Hamiltonian 

H=J2^Jb%+^-Y (2.43) 

n>0 ^ ^ 

We define the current by starting from the continuity equation. The 
current in the edge a is, by definition, given by 

dMx)=edt6p^ix) (2.44) 

and, up to a constant, we identify 

Ux) = edtMx). (2.45) 

We have then related all the quantities we want to calculate to the solutions 
of the equation of motion ()2.34|) and to the dynamics of the boson operators 
which is determined by the Hamiltonian ()2.43j) . 

2.4 Solutions of the equation of motion in 
simple cases 

It is now instructive to solve the equation of motion for the function is 
some simple cases, namely the case of two translationally invariant edges and 
when a constriction is present. 



2.4.1 Translationally invariant case 

The translationally invariant case is the simplest we can consider and we 
have for the potential the form 

Tr / l\ -,r , l\ ( VAx — X'^ Voix — X'^ \ ,^ 

= --') = [ _ J, _ ) (2.46) 

where Vi (V2) is the intra- (inter-) edge interaction^. We seek for the solution 
of the equation ()2.34j) in terms of plane waves, so we choose the form ipa{x) = 

•^Recall that the upper component corresponds to the left edge while the lower compo- 
nent to the right edge. 
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(fa{k)e . This leads us to the 2x2 eigenvalue problem 

'""[O -1 J\ipn J 2n\ V,{k) V,{k) J [ ifn J ' 

where Vi{k) and V2{k) are the Fourier transform of Vi{x) and ^^2(2^) respec- 
tively. The eigenvalues are given by 



ujk = ±\k\^Jv^{k) - Vi{k) = ±c\k\ (2.48) 
27r ' 

where c is the sound velocity. Notice that for each positive frequency there 
are two counterpropagating modes: the "up-moving" solution is 



and the "down-moving" solution is 



In these expression we have assumed A; > and we have introduced the 
edge length L. This length is assumed to be arbitrary large (as usual when 
dealing with plane-wave) and will not enter the physical results. On the 
other hand the presence of the factor l\/k is imposed by the orthonormality 
condition. The orthonormality condition also fixes u\ — v1 = l hence we can 
parameterize 

Vk = sinh(6'fe), 
where the "mixing angle" 9k is given by 

tanh2efe = :^. (2.52) 

Let us discuss briefly these results. If we have two parallel non-interacting 
edges we have 9k = for every k thus Uk = I and Vk = 0. This implies that 
the "up-moving" modes are fully concentrated on the left edge while the 
"down-moving" modes are concentrated only on the right. In this sense the 
left and right modes are well defined concepts and a good basis to discuss the 
properties of the system. Moreover we have a linear energy spectrum of these 
excitations with the edge velocity determined by the intra-edge interaction. 
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When the edges are interacting the modes propagates both on the left 
and on the right edge. The concept of left edge mode is ill-defined and 
the good basis in this case is the up-moving and down-moving modes. We 
will show, when discussing the transport, that the presence of the inter- 
edge interaction will not change the linear relation between the conductance 
and the filling fraction v when one consider the limit of low-energy (which 
corresponds to the limit /c — >^ 0) excitation. We want also point out that 
we have recovered the standard expression for the dispersion of the edge 
waves in the ordinary (non-chiral) Luttinger liquid model and that the xLL 
persists even if the interaction potential V2 is turned off. This is due to the 
anomalous commutation relation ()2.27|) between the density fluctuations on 
the same edge. 

2.4.2 A conservation law 

At this point it is interesting to discuss the existence of a conserved quantity 
for the equation of motion ()2.34|1 . We assume for the interaction potential 
the form Va,^{x — x') = Va,i3{x)6{x — x'). It is now possible to show that the 
quantity 

0l(a;X/30/3(a;) = (t>\{x)(l)L{x) - (2.53) 

is conserved 

d. (</.l(xX^</)^(x)) = 0. (2.54) 

The proof of the existence of this conservation law rests on the assumption 
of the existence of the Va^p matrix inverse"^ for every value of x. Within this 
assumption, we can consider the equation of motions for the field (pa and its 
complex conjugate 

, (2-55) 

-iuj4>\{x) = —d^(p\{x)V^,p{x)(Tl^. 

By taking the matrix Va^p to the left-hand side of these equations and then 
multiplying the first (second) equation on the left (right) by '/'jjC'I^Q, (c"a,/30«) 
and summing the results we obtain the conservation law. Notice that, be- 
cause (pa and ipa follow the same equation of motion, this conservation law 
must be verified also by ipa- 

The physical interpretation of this conservation law is interesting. It is 
possible to connect this law with the continuity equation and the conservation 

^For clear physical assumption the matrix Va.p must be symmetric. 
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of the total charge and current. We have defined the current using the 
continuity equation and up to a constant we have 



/q,(x) = —ieuj(j)a{x) 
and we can express ()2.54j) in the form 

= (lUx) - Il{x)) 



(2.56) 



(2.57) 



To appreciate the meaning of the above conservation law, let us consider the 
solution of the equation of motion. For simplicity we consider the transla- 
tionally invariant case. The left and right density fiuctuations are given by 



5pl{x) = -«e^ 

g>0 




q>0 



27rg 



2nq 

and the currents (eg = Ug 



-vbu,„ + ub^ 



t 

d,q 



-iqx 



{ubd,q - vbi g) 



(2.58) 



q>0 
q>0 




2'nq 

V 




e'''' {ub^,q + vbl^) - e-'"- {vb,,g + ublg) 
-e^'^ {ub^,q + vb\\ + e-«- {vb,^g + ublg) 



Now we define 



q>0 



q>0 V ^ 




(2.59) 
(2.60) 

(2.61) 



for the up and down moving fiuctuations densities and again from the conti- 
nuity equation 



q>0 



Ux) = -zecJ2 \h^1 V'^Kq - e-'^^&l,J 



/.(x) = -zecgy— g 




d^q 



(2.62) 
(2.63) 
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In terms of the up and down moving density fluctuations we have the simple 
relations 



c{5pu + 5pd) = c6p = lu- Id, (2.64) 

c{Spu - Spa) = I^ + h = I, (2.65) 

so that the conservation of the total charge and the total current also implies 
the conservation of the product 

{I^-h){Iu + Id). (2.66) 

Next we observe that 

5pl{x) = u5pu{x) - v5pd{x), (2.67) 

Spr{x) = uSpd{x) - vSpu{x), (2.68) 



and 



As a result we get 



h{x) \ ^ f u -v\ f Iu{x) 
Ir{x) i \-v u \ h{x) 



(2.69) 



li-iii^ It - li (2.70) 

This follows by writing the rotation matrix between [Ir, II) and (/„, Id) in 
terms of Pauli matrices and observing that 

(uao - vay)a%uao - va^) = a'. (2.71) 

Then we have reduced the conservation law written in terms of the left and 
right currents to the product of two conserved quantities for the up and 
down currents. The conservation of the currents in the up and down basis 
follows from the diagonal form of the Hamiltonian in such basis and the 
decomposition /„ — /J = (/„ — Id){Iu + Id) follows from the commutation 
rules for the up and down moving boson operators. 



2.4.3 Non-translationally invariant case 

We now consider the effect of the presence of a constriction which breaks 
the translational symmetry. This constriction can be created by depleting a 
portion of the sample by applying a voltage to the metallic gate fabricated 
on top of the mesa. When a finite k wave impinges on the constriction it can 
be partially reflected and transmitted. How this will affect the conductance 
is determined by the A: ^ limit of the reflection coefficient. If this limit is 
zero there will be no correction to the ideal Hall conductance. 
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We model the presence of the constriction by considering a piece-wise 
inter-edge potential. This choice is based on the assumptions that the inter- 
action potentials have a Coulomb origin. When the edges are forced by the 
constriction to stay close the mean distance between the density fluctuations 
is lesser than outside thus the inter-edge interaction is greater inside the con- 
striction than outside. We suppose also that the region when the inter-edge 
potential switches from the value outside the constriction to the value inside 
the constriction can be neglected. The situations we want to consider in this 
section are plotted in Fig. 12.21 where a constriction is localized in a finite 
region of the Hall bar (panel a) or the constriction extends until the drain 
(panel b). 



D 



D 



Figure 2.2: The two types of constriction we want to consider. On the left the 
constriction is localized in a finite region while on the right it is semi-infinite. 



We start by considering the case of a finite size constriction. We assume 
that the system is symmetric with respect to the a; = point and that the 
length of the constriction is d. To keep the analysis simple, we assume that 
the interaction potentials Vi and V2 are short ranged on the scale of the 
density fluctuations: this in particular implies that only points at the same 
value of the variable x interact and Va^p has the form 
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where 

( V2,i X < -d/2 
V2{x) = I V2,2 \x\ < d/2 . (2.73) 

\ ^2,3 X > d/2 

The region labelled 1, 2 and 3 are implicitely defined in the piece-wise form 
of the potential V2. We seek for a piece- wise solution. As is standard in 
scattering theory we label the full solution with the quantum numbers of the 
incident wave. An "up-moving" solution will then have the form 

f Kix)+r''^iix) X < -d/2 
^1 = 1 A-^lix) + B-^i^{x) \x\<d/2 . (2.74) 
(t^^lix) x>d/2 

The wave vectors in the regions 2 and 3 are determined in terms of the 
incident momentum ki by the conservation of the energy 

cih = C2k2 = csks (2-75) 

where Ci, C2 and C3 are the sound velocities in the three regions. In a similar 
way one can construct the solution for the "down-moving" solution 

( f^^i^ix) X < -d/2 

(pl = l A^ifl^x) + B^ipl^{x) \x\<d/2 . (2.76) 
I <(a;) +rVfc3(a;) x > d/2 

In these expressions the spinorial functions (fi^i^) fti^) given by the 
Eqs. ()2.49p and ()2.50|) where the index a has been dropped for the sake of 
notation's simplicity. 

The matching conditions are dictated by the physical requirement that 
there is not energy accumulation at the interfaces. This is equivalent to the 
requirement of the continuity of the solution at the points x = ±d/2 and 
gives four conditions to determine the coefficients A, B, t, and r for the up- 
and down- moving wave functions^. 

The solution of the set of the equations obtained by imposing the match- 
ing conditions can be obtained in a straightforward way. We get, after having 

^We must recall that the equation of motion for the spinorial wave functions ip con- 
tains only the first derivatives with respect to the time and position hence the request of 
the continuity at one point is sufficient to fully determine the solution for the scattering 
problem. 
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expressed the results in terms of the mixing angles, 



e * 2 ci 



cos{k2d) cosh(6'i — 9^) — i sm{k2d) cosh(26'2 — 6i — 63) y C3 
„ _ cos(/c2(i) sinh(^i - 63) + i sm{k2d) sinh(26'2 - 9i - 9^) .j^^^ 
cos(A;2(i) cosh(^^i — 9^) — i sm{k2d) cosh(26'2 — 9i — 9^) 

j^u cosh(6'3 - 6*2)6-*^^°' /c7 



(2.77) 



5" 



cos{k2d) cosh(0i — 9^) — i sin(fc2(i) cosh(2^^2 — 9i — 9-3) y C2 ' 
sinh(^3 - ^2)e*^'^ /cT 



cos(A;2(i) cosh(^i — 9^) — i sin(A;2(i) cosh(26'2 — 6'i — 6'3) y C2 



The coefficients for the down-moving solution can be obtained by the substi- 
tutions 

t« ^ t'^- r" ^ r"^; A" B'^; A'^, (2.78) 

and 

9i ^ ^3; ^^3 -> ^^i; ci ^ C3. (2.79) 



It is easy to verify that the wave function p.49|l with the coefficients deter- 
mined by ()2.77|1 satisfies the conservation law ()2.54|) . 

We recover the case of a semi-infinite constriction by considering the limit 
c? — >■ in the expressions for the transmission and reflection coefficients 



cosh(^i -^3)Y C2' , . 

„ ^ sinh(g, - ^3) ^'-'^^ 
^ cosh(^i - ^3) ■ 

If one defines the interaction renormalized filling factor Ui = z/e~^^' for the 
various region then it is possible to rewrite the above results as 

^1 + ^3 ' , , 

o r- 2.81 



^i + ^sy C3 

These expressions remain valid when we consider the case of two regions with 
different filling factors [l^. In the following we will consider the symmetric 
case 9i = 9^ i.e. the interaction is symmetric with respect to the center of 
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the constriction. From the expressions ()2.77|) we get 



^—ikd 



cos{k2d) — i sm{k2d) cosh(26'2 — 26 



2sin(M)sinh(2e2-20i) .^^^ 

6 1 



fci+fc2 ^ I — (2.82) 



A"" 



cos{k2d) — i sin(/c2'^) cosh(26'2 — 26i 

cosh(gi - 62)6-'^'^ fcl 
cos{k2d) — i sm{k2d) cosh(26'2 — 26*1) y C2 ' 

sinh(ei - 62)6''^'^ /cT 



cos{k2d) — i sm{k2d) cosh(26'2 — 26*1) y C2 



We will use these expressions to calculate the Hall conductance in the pres- 
ence of a constriction. 

The generalizations to the case when many constrictions are present or a 
constriction connects regions with different filling factors are straightforward. 



2.5 The tunneling Hamiltonian 

We want now to consider the effect of the tunneling between different edges. 
The physical origin of tunneling lies in the fact that the electron quasi- 
particles are not completely localized in one or the other edge i.e. the density 
matrix p(r, r') has a finite value even if r 7^ r'. 

The physics of the tunneling is obviously lost in the hydrodynamical ap- 
proximation. We need then to insert the tunneling by hand in our Hamilto- 
nian. To do that we need to define a quasi-particle creation operator which 
adds a quasi-particle with charge e* (not necessarily equal to the electron 
charge e) at the point x of the edge a. This is accomplished by requiring 
that this operator "^^ix) satisfies the commutation relation with the quasi- 
particle density 

i^iix), Sp^ix')] = -^So.,pS{x - x')^i{x). (2.83) 

At the best of our knowledge there is not exist a general theory which predict 
the correct value of e* for arbitrary value of the filling factor. For certain 
values of the filling factor, as an example those given hj u = l/(2m+ 1) where 
m is an integer, it is believed that e* = ve. Two approaches are then possible. 
One can fix e* = ve and then makes a comparison with the experimental 
results. On the other hand it is possible to treat e* as a phenomenological 
parameter determined by the confrontation with the experiment. 
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The equation ()2.83|) allows us to express the quasi-particle creation oper- 
ator in terms of the boson operator via 



exp 



X exp 



n>0 
n>0 



V 



(2.84) 



where the unitary fermion operator^ f/J^ commutes with all the boson oper- 
ators and increases the total charge by — e* 



where the total charge operator is defined as 



(2.85) 



Qc 



(2.86) 



The solution of the equation ()2.83|) can be obtained by the observation that 
the commutation relation can be converted to a differential equation by using 
the different representations 



which in turn imply that 



d_ 



dg{b) 
dh 



(2.87) 



(2.88) 



The operator Ua is then the "arbitrary constant" in the solution of the differ- 
ential equation and its properties can be deduced from the physical insight 
that it must increases the total charge Qa- We have also written the quasi- 
particle creation operator in a normal ordered way: this will avoid some 
complications with the normalization [s^ . 



^We have omitted a normalization constant which depends on a short-range cut-ofT. 
We choose this normahzation constant such that this operator is dimensionless. 
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In terms of the quasi-particle operators, the tunnehng between the edges 
coupled by a constriction at x = is described by the Hamiltonian 

Ht^T: *[(0)*r(0) : +T* : ^^,(0)^^(0) : (2.89) 

where F is the (phenomenological) tunnehng amphtude and the symbol : . . . : 
indicates the normal ordering. The complete Hamiltonian then reads 

H = Ho + Ht = Y, ^nh% + r : ^i(0)^K(0) : +r* : ^1,(0)^^(0) : (2.90) 

n>0 

where now the tunneling Hamiltonian can be viewed as an interaction Hamil- 
tonian for the bosons and now the total charge in a given edge is not a con- 
stant of motion: its time derivative defines a tunnel current operator It as 
follows 

(2.91) 

= f-{V : *U0)*«(0) : -r : ^Uo)^l(O) :). 

In the following chapter, when wc will deal with the formulation of the 
transport and derive the relations between the current and voltage, we will 
describe a method to calculate this tunneling current by defining an exact 
boson propagator and developing a perturbative scheme to evaluate its ex- 
pectation value. 

We must point out also that we do not specify any statistics for the quasi- 
particle operator. Wen was able to show that the quasi-particle operator 
follows a Fermi algebra if and only if v is the inverse of an odd integer. If 
V is the inverse of an even integer the quasi-particle follows a Bose algebra. 
In the intermediate case the quasi-particle has not a definite statistics. We 
do not restrict the range of variation of the filling fraction and then we do 
not have real fermion operator for the quasi-particle. We will show in the 
following that the restriction to fermion operator is not necessary and we can 
fully develop a perturbative approach to the tunneling. 

2.6 The multi-probe setup 

In the previous sections of this chapter we have introduced our model for the 
quantum Hall bar and the way to describe the edges and their excitations. 

For the sake of simplicity in introducing the notation and the various con- 
cepts we have treated the case of two parallel edges with only two probes. 
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namely the drain and the source. However as can be seen from the scheme 
of the experimental setup in Fig. ll.3l we need to consider a more complicated 
situation where many probes, edges and contacts are present. To do that we 
consider the schematic in the Fig. 12. HI where we have introduced six different 
probes, two of them are used to inject a steady-state current in the system 
(the D and S probes in the figure) and the others are used to measure either 
the Hall or the longitudinal voltage. The presence of the constriction is also 
included. 
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Figure 2.3: A schematization of the experimental device where many edges 
and probes are represented. The gates (G) create a depletion zone forbidden 
for the electrons and then force the edges to stay close. 



We choose the following notation. Every edge is identified by an arable 
figure, while the probes are indicated by roman figures. The point where 
an edge exits the device is determined by the notation Xa,m where a is the 
arable figure which refers to the edge, while m identifies the probe. Generally 
a point belonging to a given edge is identified by its position Xa with respect 
to the edge. The number Xa varies in [—00, 00] but we consider that when 
the electron leaves the device the correlation vanishes exponentially with the 
increase of the position. 

We need to change also the definition of the density and quasi-particle 
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operators. The density fluctuation that belongs to the edge a is now indicated 
by 5p{xa) while the quasi-particle creation operator is ^^(xq). The simplest 
way to think at this change of notation is to attach the index a directly to the 
variable x rather than to the operator. The usefulness of such notation will 
be proved when we will deal, in the next chapter, with the calculation of the 
transport properties. When we consider the delta function of the position, 
namely the quantity 5{xa — xp), we identify it with the product of a delta 
function relative to the edge index times the delta function of the position 

5{x^-xp) = 5^,p5{x-x'). (2.92) 

The intra-edge interaction is substantially not affected by this change 
of notation. For the inter-edge interaction we consider that our device is 
specular with respect to the In such a way when we consider the 

interaction potential V2{xa—x'p) we consider two point with the same abscissa 
belonging to two different edges. 
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Chapter 3 



The formulation of the 
transport 



In this chapter we derive the expression for the conductance matrix in terms 
of the correlation function of the quasi-particles. In the transport theory we 
need to calculate the current induced into the reservoir i due to a change 
in the potential Vj of the reservoir j. If we consider only the linear relation 
between the change in the potential 5Vj and the induced current 5Ii we have 



As it is recalled in this definition, the conductance matrix Gij is a function of 
the potential of all the reservoirs. In the following, to simplify the notation, 
this dependence will be understood even if not indicated explicitely. We 
choose the convention that a current is defined positive when it enters a 
reservoir (leaving the system) and negative viceversa. The matrix element 
Gij are also subjected to the physical conditions 



which represent the gauge invariance and the charge conservation, respec- 
tively. These expressions determine the value of the diagonal elements On 
once the off-diagonal elements are known. Notice also that in general there is 
not other way to calculate these diagonal elements. The rest of the thesis will 
deal with the calculation of the conductance matrix in various cases study. 




(3.1) 




(3.2) 
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3.1 The conductance matrix 



We want now to express the conductance matrix in terms of the boson cor- 
relation function. 

The starting point is the definition of the current of the edge a in terms 
of the displacement field 0. As we have discussed in the previous chapter 
this relation can be derived from the continuity equation which gives 

I{xo,) = edt(p{xa). (3.3) 

The current in the terminal i is then the algebraic summation of all the 
currents that enter or leave this terminal. We define the function 

{+1 if q; enters i 
—1 if a leaves i (3.4) 
otherwise 

such that the current on the terminal i is given by 

h^^^a,iH^a,i)- (3.5) 
a 

This relation and the function ^ represent the mathematical formulation of 
our convention on the sign of the currents flowing in the system. In general 
this is the sum of the two different currents, one leaving the terminal and 
the other entering it. We choose, as the versus of the total current, the 
versus from the source to the drain in the two edges. This choice fixes the 
"contact function" j and, as we will see, the sign of the potential drops at 
the contacts. 

By remembering that in the linear response theory the voltage is coupled 
to the charge density, we have 

oo oo 

SI{x^,uj) = tjYl J dx'f,5V{x'f,) J dt{[dt<P{xa,t),d,^,(j){x0')])e'^^+'^^' 

^ -oo 

(3.6) 

where uj is the frequency of the external potential and (. . .) is the equilib- 
rium average. We are interested in the static limit hence in the following 
expressions the limit cu ^ 0, when not explicitely indicated, is understood. 
After an integration by parts with respect to the time, we arrive at 

5I(xc,) / dx'f^ [<^9^'i3D{xa,x'f^;u;)]5V(x'f)), (3.7) 
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where 



oo 



2m /([0(x,,t),0(x')])e^(-+^'^)*c?t 



(3.8) 





is the retarded displacement-field propagator, whose explicit expression in 
terms of "phonon eigenfunctions" is 



d^^>D{xa, Xfs'; t = 0+) = -27r[0(x„), 9^^/0(x/)] oc 5(x„ - x^') (3.10) 



vanishes unless Xa = xp'. We will show in the following that the evaluation 
of this function on the point Xa^i is related to the conductance. However, 
as we have discussed briefiy before, the conductance matrix elements we can 
calculate are the off-diagonal ones only. Hence in doing these calculations 
we deal only with the 7^ xp' and this validates our procedure of 

integrating by parts. 

We need now a model for the external potentials which are applied to the 
reservoirs. We model only the potentials in the region outside the device. To 
do that we consider that the potential is uniform in the portion of the edge 
that runs inside the reservoirs, drops to zero at the contact points between 
the leads and the device, and is zero inside the device. Thus we choose an 
external potential which satisfies 



where the "contact function" was defined in the Eq. ()3.4j) . This expression for 
the external potential makes very easy to calculate the conductance matrix. 
Indeed after an integration by parts in the Eq. ()3.7|) . carried out on the 
assumption that the correlation function decay exponentially for x ±oo, 
we obtain the expression of the current arriving at the reservoir i via the 
edge channel a due to a potential variation applied to the edge channel j3 by 
the reservoir j {i ^ j) 




(3.9) 



In performing the integration by parts we have used the fact that 





j 




(3.12) 
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We have thus obtained an expression for the conductance matrix in terms of 
the boson correlation function. 

As an apphcation we want to calculate this function in the case of two 
translationally invariant edges. In this case the relation between the energy 
and the momentum is uj = kc where c is the phonon velocity and the wave 
functions are given by the expressions ()2.49j) and ()2.5()j) . To ensure that the 
wave function vanishes in the limit |x| ^ oo we shift, in the argument of the 
complex exponential, k ^ k + irisign{x) in the ()2.49|) and k ^ k — irjsignlx) 
in the ()2.50|) . After plugging in these wavefunctions in the expression for the 
correlation function we reduce the sum over k > to the integral over the 
real axis in k. We readily see that only the poles ck = ±{uj + irj) contribute 
to the integral thus giving 



-i limujD^XayXp' ,t) =i'Q{xa — x'^) 

+ uQi^x'p — X 



(3.13) 



—UV 
,2 



f ^ —UV 
—UV u^ 



where the functions u and v are the limits — >^ of the functions Uk and Vk 
we have defined in the ()2.51|) . In the case of non-interacting edges we have 
u = 1 and u = 0, thus we have upward propagation on the left edge and 
downward propagation on the right edge. This makes the conductance Gij 
vanishes unless the reservoir j is not "upstream" the reservoir i according to 
the definition of an ideal Quantum Hall system. 

When the edges are interacting, from the direct application of the Eqs. 
(jTT^ and dSini), we obtain 

G21 = G12 = (3.14) 

Notice that in this simple case an edge contacts both the reservoirs. One 
might wonder how this result modifies the relation between the current and 
the source and drain potentials. As pointed out by Wen |lj], when the 
interaction is present, the edge potentials are a linear combination of the 
source and drain potentials with the coefficients determined by the mixing 
angle 9. By taking into account this relation we arrive at 

I = j,''1{Vs-Vd) (3.15) 

and this implies that a measure of the Quantum Hall conductance is insen- 
sitive to the presence of a translationally invariant inter-edge interaction. 
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A similar calculation can be developed for the case of a constriction by 
using the reflection and transmission coefficients in Eq. ()2.77|) . In this case 
the calculation is complicated by the piece-wise form of the wave functions. 
To keep the analysis as simple as possible we consider the case when four 
reservoirs are attached just above and below the constriction and the interac- 
tion is present only inside the constriction region (i.e. we fix 6'i = ^3 = 0). By 
moving from the wave functions ()2.49|) and ()2.5U|) we have obtained fixing, 
as an example, the first reservoir 



00 



g2 !• Qik{xi-X2) 



-00 



Gsi = 0, (3.16) 

\ r ^-ik{xi+X4) 

G41 = u- / dk- ; -r^ 

h27ri J k- u/c - iO+ ^ 



With ^1 = ^3 = the reflection and transmission coefficients read 

'^k 



cos{k2d) - i sm{k2d) cosh(26'2) ' 

ism{k2d)smh{292) ^.j^.d 
cos{k2d) — zsin(A;2'^) cosh(2^^2 



-e 



where we have Ciki = 02/^2- In performing the integration to calculate the 
conductances we observe that we must close both the integrals in the upper 
half complex plane. Since tk and have only poles inside the integration 
path in the lower half complex plane, the only pole is at k = uj/c + irj and 
we obtain 

G21 = i^^t^iuj/c), 

\ (3.18) 

h 

The ideal Quantum Hall conductance is recovered by observing that in the 
limit of vanishing frequency ^ we have — 1 and — > 0. We are 
then arrived to the conclusion that the constriction and the interaction do 
not modify the ideal Quantum Hall conductance. Indeed one must observe 
that the interaction does not allow for the passage of one quasi-particle from 
one edge to the other. On the other hand the constriction, in the limit of 
small energy, becomes fully transparent thus recovering the translationally 
invariant results. 
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3.2 The tunneling conductance 



When we consider the tunnehng, we know that the full Hamiltonian cannot 
be diagonalizable. Thus we develop a perturbative scheme to calculate the 
tunneling conductance i.e. the correction to the ideal Quantum Hall conduc- 
tance induced by the presence of the tunneling. We choose to develop this 
theory starting from the equations of motion. Such an approach does not 
rest on any assumptions about the statistics of the quasi-particles and, then, 
can be applied to the case when v is not the inverse of an odd number. 

As we have pointed out the presence of the tunneling adds to the free 
Hamiltonian an interaction potential between the bosons operators. The task 
at hand is then to calculate the corrections, due to this interaction, to the 
displacement field 0. We show that this correction can be exactly expressed 
in terms of a tunneling current propagator and we provide a perturbative 
evaluation of the latter. 

We introduce now some compact notation. We define the phonon spino- 
rial operator 

K ^ ( b ) (3.19) 



where the index i G {1,2} (the upper part corresponds to 1), and the asso- 
ciated phonon propagator 

D'rfAt)^-l^(t)i[K{t),Bi]]). (3.20) 
We define also the spinorial eigenfunction 

^ ( Xtr ) • P-^^) 

With these definitions the phonon field propagator can be written as 

i,j {n,n'}>0 

In the following we will suppress the summation symbols and the sums over 
repeated indices will be understood. The phonon operator satisfies the equa- 
tion of motion 

id,Bi, = nifBi - ^It (3.23) 

where we have defined 

^if-("'n ° V ^: = f V (3-24) 
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and 

1: 



^- = \j—Y.^n{^-)- (3.25) 



It is now straightforward to verify that the phonon propagator satisfies 
the equation of motion 

{idA,i - ^M]^^ = ^-^kA,nm - ^QUt) (3.26) 
where the auxihary operator 

Q'nit) = -'-m{[lT{t\Bl}]) (3.27) 

{idA,, - = -^MAt) (3.28) 

*2 

Mrit) = Mrit) - ^{HT)S{t) (3.29) 

MT{t) = -^&{t){[lT{t),lT]) (3.30) 

is the tunnehng current propagator. We will show later that this propagator 
can be expressed in terms of the quasi-particle correlation functions. 



in turn satisfies 



with 



where 



The Eqs. ()3.26|1 and ()3.28|) can be readily solved by Fourier transforma- 
tion to obtain 



n.n' L sn.n' 



* (3-31) 

+ Mt(c.)(f™)* {[d^'^]:;^m) 

where [D^^^lf^, (u) is the noninteracting phonon propagator. Because Mt is 
related to the quasi-particle correlation functions, we have then expressed the 
phonon propagator in terms of the tunneling current propagator. The non- 
interacting phonon propagator is determined by the Fourier transformation 
of the operator 

[D^"^]if,n'i^) = ^-^kA^mW^, - ^)-\ (3.32) 

Notice that there are no approximations on the calculation of the tunneling 
amplitude. Thus this set of equations constitutes an exact approach to the 
problem of the tunneling. 



55 



Let us denote by Cf^ the conductance obtained in the absence of tunnel- 



ing and by 



5Gij — Gij — (3.33) 

the correction due to the tunnehng. By starting from the definition ()3.12|) . 
with some straightforward manipulations, one arrives at 



^ hm > U[^'"1la;m,U^;^) 

(3.34) 

X [a;MT(^)]5^[Z^(°t(0.,4.;^)^/3.- 

where the indices 7 and 8 run over the two edges that are coupled by tunneling 
at a; = 0, and the Green's function of the noninteracting displacement field, 
afi[x-,x' \uj\ is given by Eq. ()3.9|) . 
As a concrete example, let us consider a four-terminal geometry, as may 
be obtained from Fig. 12.31 by considering only the terminals from 1 to 4. Let 
us assume for simplicity that the mixing angle Q is independent of x. Then 
from Eq. ()3.13|) we immediately get 



where the upper (lower) component refers to the left (right) edge and u = 
cosh 6', V = sinh6'. 

Substituting this in Eq. ()3.34p we find 

"^^^i = ^Go^^^^J i^i^ (3.36) 



where 



Mt{uj) f , / —uv 



6Gap{x,x') = - ie-^^ lim ' { e(x)0(-x' 



,2 



f ^ —UV \ ^ / X ^ / /X / —UV U^ 



+ Q{-x)Q{x') 2 + Q{x)Q{x 



,2 



—UV U I ^ ' ^ ' \ V —UV 

+ e(-x)0(-x') ■ ^ ^ 



—UV V 



U^ —UV 



(3.37) 
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Putting this in Eq. ()3.36|) and noting that ^q,i = = 1, = ^/33 = — 1 
(with the labels i, j as specified in the figure) we finally obtain the correction 
to the Landauer-Biittiker conductances of the ideal system: 



5G. 



le 



hni 



uv 



V 



—uv —V 

2 



2 

UV 



—UV 



UV 



uv —u —uv 

2 



V 
UV 



(3.38) 



In the following for simplicity of notation we define the function 



9t 



i lim 



Mt{uj) 



UJ 



(3.39) 



We consider the specific setup of Fig. 12.31 which represent an ideal rep- 
resentation of the experimental setup proposed by Roddaro et al. [is!]. The 
resistance Rxx of the quantum point contact is measured between terminals 
3 and 4 



(3.40) 



where I is the source-to-drain current. By considering that the constriction 
does not affect the source and drain probes, the full conductance matrix reads 





1 














-1 


\ 







1 





-1 












-1 





1 + 5gn 


^9l2 


S913 


S9u 












-1 + 5g2i 


1 + ^922 


S923 


5 924: 









-1 


^931 


5 932 


1 + 5g33 


S934 




v 








^9ii 


S942 


-1 + Sg43 


1 + 55(44 


/ 



(3.41) 

where the indices i, j run over {S,D, 1,2,3,4} and the right bottom sub- 
matrix is given by 6gij = In the next chapter we will evaluate the 
gT in a perturbative way where the small parameter is the tunneling am- 
plitude V: we will show that the leading term in 5Gij is proportional to 
|rp hence it must constitute a small correction to the ideal Quantum Hall 
conductance. In the Appendix [0 we discuss briefly the Landauer-Biittiker 
formalism and give some insights about the calculation of this matrix. As it 
is customary in the experimental setup we flx Vd = 0, Is = —Id = —I and 
Ii = I2 = I3 = h = 0. With these constraints, the equation 1)3.11) can be 
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easily solved^, and to the lowest non vanishing order in we get 



Rxx = -^e-^\u + v)u qt = -^e~^ cosh 9 qt- (3.42) 

Notice that this perturbative result is valid only so long as R^x is much 
smaller that ^: The tunneling amplitude V must be sufficiently small to 
justify this perturbative approach. 



3.3 The tunneling amplitude 

In the previous section we have connected the conductance matrix to the 
current-current correlation function Mt- The next step is to connect this 
function with the quasi-particle correlation function. In this step we intro- 
duce a perturbative approach and develop all the terms to the second order 
in the tunneling amplitude F. 

In Mt two terms are present. The first one is the current-current correla- 
tion function proportional to {[Irit) , It]) ■ Because the tunneling current It 
is proportional to |r| this correlation function is already second order in the 
tunneling amplitude. The second term is the average on the ground state of 
the operator Ht- This average must be calculated to the second order in |r|. 
By using the Kubo formula ^22] we arrive at 

t 

{HT{t))='- I dt' {[HT{t'),HTm- (3.43) 

— oo 

To simplify the notation we define the operator 

A{t)=:^[{0,t)^n{0,t): (3.44) 
and then write the tunneling Hamiltonian as 

HT{t) = rA{t) + r*A\t) (3.45) 
and the tunneling current as 

/T(t) = ^(rA(t)-rMt(t)) (3.46) 

"'^ Notice that the gauge invariance and the charge conservation constraints make singular 
this matrix. We need then to ehminate a row and a column before proceed with the 
standard technique to solve the linear system. 
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where we have used the fact that the tunnehng takes place at the point x = 0. 
In this notation we have then 



{[HritrHTm = \rmA{t),A\t')]) + {[A\t),A{t')])) (3.47) 

and 

{imjAo)]) = '-^\mm),A'm) + {[A\t),Am)). (3.48) 

Notice that in doing these averages the anomalous means that involved two 
A or two operator have been dropped. This is justified by the presence 
of the operator Ua in the definition of the quasi-particle operator. Indeed, 
being unitary operators, we have U\^Ur = U\Ul = 1 while the averages 
of the others product of two U operators is zero. We have restored here 
the left L and right it! notation for the sake of simplicity. Indeed in the 
constriction only two edges are present (we do not have any contact inside 
the constriction). Notice also that we have chosen the coordinate systems in 
these edges to have the a; = point at the same height of the sample. 

The simple relation between a commutator and its hermitian conjugate 

[A,B]^ = -[A\B% (3.49) 

if we define the correlation functions 

G^{t-t')^{A{t')A\t)), 

G^{t-t') = {A\t)A{t')), ^ ■ ' 

allows us to recast the average of the tunneling Hamiltonians commutator as 



{[HT{t%HTm =|rp {{[A\t),A{t')]) - {[A\t),A{t'W) 

=2i\T\Hm{G-{t']t)-G+{t']t)) (3.51) 
=Ai\V\HmG_{t']t) 

where we have used the relation G^{t'\t) = G*_{t';t) that we will prove in 
the following chapter. The tunneling current correlation function can be 
rewritten as 

*2 

([JT(t),/r]) = -4z^|r|2lmG'_(t;0). (3.52) 
Because by definition we have 

Mt{u)^ Mt{u)-^ f dt 5{t){HT{t))e'''^ (3.53) 
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we must calculate the Fourier transform of the tunneling correlation function 
and of the second term in the right hand side. To calculate this second term 
we use the Dirac S to perform the integration on time and we have 



oo 

/ 



dt 5(t){HT{t))e"^^ = -4 



4- 


h 


4- 






h 




c 

' 1 



— oo 

oo 



h 



dt' ImG_(i';0) 
W ImG'_(-i';0 
dt' lmG^{t';Q) 



(3.54) 



where we have used the fact G-{t;0) = G^{t;0) = G+(— 1;0) we will prove 
in the next chapter. It is now easy to obtain the relation between Mt{uj) 
and the function G-{t). We get 



lim 



Mt{oj) 



4e*^|r 



oo 

lim / 



dt 



cos{ujt) 



'-luiG_{t) 



lim 



oo 

/ 



dt — ^ImG_(t) 



(3.55) 



4ie*^|r 



dt tlmG_{t). 



The presence of a bias voltage on the edge a modifies the time evolu- 
tion of the corresponding quasi-particle operator 



(3.56) 



The underlying physical assumption is, of course, that each edge is in equi- 
librium with a reservoir, from which it generates, at potential V^. Under 
this assumption, the bias voltage dependence of the conductances can be 
calculated with no additional effort. Indeed the function can be easily 
expressed in terms of the quasi-particle creation and annihilation operators. 
Prom the definition we have 



G.{t-t') = {A{t')A\t)) 

= (: *U0,O*ii(0,O :: ^^^^ lM :) 



(3.57) 
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thus this correlation function is modified by a phase factor due to the presence 
of the finite edge voltages. Notice that the same phase factor modifies the 
function G+. 

With arguments similar to those used to derive the Eq. ()H.55j) we have 

{[HT{t'),HT{t)]) = 4icos{ujTt')lmG4t';t), 

4ze*2 , ^ , ^ (3.58) 



{[lT{t),lT]) = ^cos{ujTt)lmG40;t), 



and 



2 2°° 

(yfyfcuy) = i lim — = — [ ^ [ dt t cos(ujTt)lmG ^(t; 0) 
i^^o uj n-^ J 

oo 

j dt sin(u;rt)ImG_(t;0) (3.59) 





oo 

4e*2|rP d 



dujT 



oo 

4e*2|rP d 



where we have defined the frequency 



Im j dt e*'^^*ImG'_(t;0), 



ujt = — —. (3.60) 
n 

We evaluate the tunneling voltage, the voltage across the quantum point 
contact, by 

Vt = + (3.61) 

and it is possible to prove that Vr = Vh where Vh is the Hall voltage. 
The amplitude gTi^Jr) niust be plugged in the Eq. ()3.38p to calculate the 
correction to the ideal conductance. The result of this operation is the main 
topic of the next chapter. 
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Chapter 4 



The transport properties in the 
presence of the constriction 



In this chapter we apply the formahsm we have developed in the previous 
chapter to the case when the tunneling takes place inside the constriction. 
We need to calculate the retarded response function G-{xa, xp', t) defined 

as 

G_{x,t;x',t') = ¥^{x',t')<iln{x',t') :: ^^x, t) :) . (4.1) 

To do this calculation we use the Haussdorf lemma 

= (4.2) 

and the expression of the operator in terms of the boson operators as 
defined in Eq. 

We start by considering the zero temperature limit and obtain an approx- 
imated form of this correlation function. The presence of the constriction in- 
troduces the time scale, to = d/c, i.e. the time a density fluctuation needs to 
travel through the constriction. We can then write the correlation function 
as the sum of two different contributions, one for t < to and the other for 
t > to- This allows us to perform the Fourier transform we need to calculate 
the current. 

The non-zero temperature case is then recovered by a conformal trans- 
formation [3^ that allows us to obtain the form of the correlation function 
directly from the zero temperature expression. Also in this case the corre- 
lation function can be separated in two different contributions for short and 
long times. 
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4.1 General results 



By plugging in the definition of the G-{x,t; x' ,t') correlation function the 
general form of the quasi-particle operator \E'^(xq,) given by the Eq. ()2.84|1 

we get 



G-{x, t] x' , t') = exp 



27re 



*2 



n>0 



(4.3) 



x««(x') + ^:x(a:'))e^^'=(*-*') 
while when we perform the same calculation on G^{x, t; x', t') we have 

27re*2 



G+(x, t; x', t') = exp 



n>0 



(4.4) 



To obtain these relations we have used the fact that the temperature is zero. 
In this case the average of the boson operator exponentials on the ground 
state gives 

^gafetg/36„^ = 1 (4.5) 

where a and (3 are arbitrary coefficients. We will see that in the case of 
non-zero temperature this average is a function of the temperature via the 
Bose distribution. 

From these expressions we get immediately the relations 



G4x,t;x\t') = Gl{x,t;x',t') 
and if we fix a; = = we have also 

G-{0,t;0,0) = G_(t) = G+{-t). 



(4.6) 



(4.7) 



This validates our results obtained in the previous chapter about the relation 
between the tunneling amplitude and the G- {t) correlation function. Notice 
that now and in the following we use the fact that the times t and t' appear 
always in the combination t — t' to set t' = 0. When we fix a; = = we 
obtain 



G- (t) = exp 



27ie* 



2 iuj^t 



n>0 



(4. 
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To calculate the effects of the constriction we substitute in the Eq. ()4.3p 
the wave functions ip determined in this case. The result of this operation is 



exp 



fc2>0 



(4.9) 



u —ik'zx 1 2 



d —ik2X 



d ik2X 1 2 



where the coefficients A", B^, A'^ and 5'^ are given by ()2.77p . By assuming 
that the tunneling is localized at the point x = x' = Q and that the system 
is symmetric with respect to the center of the constriction, i.e. we assume 
Oi = 03, and substituting the expression ()2.77|1 in this equation we obtain 
our key result 



G^{t) =exp 



Ane*^ cosh(2g2) ^262 
Lue^ cosh(2ei)^ 



E 

k2>0 



^ik2C2t 



X 



cosh(2^i2) — sinh(26'i2) cos(A;2(i) 



1 + 2sinh2(^i2)sin2(A;2rf) 
where we have defined 612 = 61 — 62. 



(4.10) 



4.2 The translationally invariant system 



Let us discuss, starting from the Eq. ()4.10j) the case of translationally invari- 
ant edges. To do that we fix 61 = 62 = 6. With this position the correlation 
function ()4.10|1 greatly simplifies to give 



G_(t) =exp 



Aire 



* 2 



-26 



Lve^ 



fc>0 



Akct 



(4.11) 



To calculate this expression we use the well known analytical results 



E 

n=l 

00 

E 

n=l 



cos(ng) 
n 

sin(ng) 



n 



which can be summarized in 



E 



n=l 



n 



iln(2-2cos(g)) 



- In (1 - e^^) . 



(4.12) 



(4.13) 
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If q is small compared to 1 we have the approximated results 



^^^-ln(g), E^r^^2' E— (4-14) 

71=1 n=l 71=1 

To calculate the expression ()4.1H) we write the momentum k in terms of a 
integer j defined hy k = 2nj/L and then evaluate the series in to obtain 



G- (t) = exp 



2e* ^29 



■e 



/ Indict r\ 

In M - e—-^\ 



(4.15) 



and in the limit of large system size ct/L '^1 we have 



G_W=(«-^) (4.16) 

where 5 assures the convergence of the series even when t ^ 0. This function 
is the propagator for the Luttinger Liquid model, with the anomalous expo- 
nent -^e~^^. Notice that if we assume e* = ue we get for this exponent 



2z/e-^'' = 2z>. (4.17) 

In this case the tunneling differential conductance at zero temperature is 
predicted to have a power law behavior with exponent given by 2(e*/ e)^/ i/— 2. 
In fact to obtain the relation between the tunneling current and the tunneling 
voltage we must evaluate the Eq. ()H.59j) . To perform this calculation we 
introduce a Heaviside function, 0(t), to extend the integral over the real axis 
and then we evaluate the Fourier transform by calculating the convolution 
between the Fourier transforms of the G±{t) and of the 0(t). Some details 
on the calculation of Fourier transform of the G±{t) are discussed in the 
Appendix O The result of these operations is the tunneling amplitude 

( ^ / ^ V ^ I 1-1 ■ i ^ 



where we have defined 



2 P*2 

a = - e-'. (4.19) 



The Eq. (I4.18|l reproduces the result of Wen for the tunneling at zero tem- 
perature in the translationally invariant case 0] . We want to point out that 
in the calculation of Wen a ultraviolet cut-off determined by the magnetic 



^Notice that the lower hmit fc > forces the sum on j to start from 1. 
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length scale appears. This cut-off is necessary if one evaluates the sum in 
the correlation function as an integral to guarantee the convergence. In our 
case however an infrared cut-off determined by the sample length L appears. 
The finiteness of this length assures the convergence of the series in the G- 
correlation function eliminating the singularity at A; — 0. 

When we consider the finite temperature case we must evaluate the boson 
thermal average {h\Jjn) which appears in the argument of the exponential in 
the definition of the G- correlation function. The result of this evaluation 
is that the correlation function at finite temperature can be obtained from 
that at zero temperature by the conformal transformation [ssl 

(5±,t)^^^^^i!^^|±M. (4.20) 

Notice that we are using units in which h = ks = 1, where ks is the Boltz- 
mann constant. The correct physical dimensions are restored via the substi- 
tution T — > ksT/h and this is understood in the following equations. With 
this substitution the evaluation of the qt function follows the same steps we 
have discussed in the zero temperature case. In this case however one must 
evaluate the integral 

j„, (■'TtT 



The evaluation of this integral is discussed in the Appendix O We have thus 
obtained 



gT[ujT,T) =2^ — — — - — (ttT) sm ' — i 



\2ticJ ^ ' V 2 / 

d f a ujt 

miB I — — , 1 — a 



(4.22) 



Out V2 27rT' 
where B is the Euler Beta function defined by [l^l 

,4.23, 

With some straightforward manipulations on the B function we can recover 



the result of Wen [1J| for the finite temperature tunneling conductance. 



4.3 The effect of the constriction 

We consider now the effect of the constriction which breaks the translational 
invariance. To do that we must evaluate the sum which compares in the G- 
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correlation function in the Eq. ()4.10|) . We define the quantity 
27rcosh(2^2) 



S4t,d) 



L cosh(2^i) 

eifc2C2t /cosh(20i2) - sinh(20i2) cos{k2d) 



X 



E 



^^^P h V l + 2sinh'(^i2)sin2(M) 



(4.24) 



and in the following we will evaluate it numerically. We then separates its 
real and imaginary parts 

jl^g _ ^ cosh(26'i2) - sinh(26'i2) cos{2TTdn/ L2) 



X 



1 + 2 sinh^ 6112 sin^(27rc/n/L2 
cos(27rnc2t/Iv2) 



ImS_{t,d) =J2 



n 

cosh(26'i2) — sinh(20i2) cos{2TTdn/ L2) 



(4.25) 



1 + 2 sinh^ 012 sin^ {2tt dn / L2) 

sin{2'KC2tn/ L2) 
n 

Notice that in these functions we have defined the integer n such that k2 = 
2Tcn/L2 and the length L2, defined by 

^ cosh(gi) 

L2 = L — — — , 4.26 
cosh(6'2) 

which takes into account the different propagation velocities of the wave 
function in the regions 1 and 2. The convergence of the series is guaranteed 
by the oscillatory behavior of the trigonometric functions. 

The constriction introduces the characteristic time to = d/c2 which is the 
time an edge wave needs to travel through the constriction. We discuss the 
two limits of t ^ to aiid t <^ to- Iii both these two limits we recover an ex- 
pression for the G- correlation function which is similar to the translationally 
invariant case and this allows us to obtain the tunneling conductance. 

First let us consider the c? — regime. From the Eq. ()4.25|1 we get 

ReS^it.d^Q) =e-2^i2 



2 

e-2^iMn(27rc2t/L2), (4.27) 



ImS.{t, ^ 0) =e-^^i2 {^{l - Ac2t/L2 
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By plugging in these results in the expression for the correlation function 
we get 



it) 



6 



2TTiC2t 



From the definition ()4.2(j|l and from the energy conservation kiCi 
have also the relation 

L2 L 



(4.28) 
/C2C2 we 

(4.29) 

C2 Cl 

With this result we recover the same correlation function of the translation- 
ally invariant case when the constriction is not present. 

In the other limit (i — 00 we have substituted in the two functions in the 
Eq. ()4.25|) the averaged values, (cos(/c2'^)) = 0, (sin^(fc2'^)) = 1/2 obtaining 

ReS-{t,d^ 00) =(1 + tanh^(^i2)) 



X 



In 2 - 2 cos 



/ 27rc2t\ 



c:^- {1 + tanh^ 
ImS^{t, d — > 00) =(1 + tanh^ 

~(1 + tanh^ 



\ L2 J 
12)) ln(27rc2t/i:2) 

2C2t/L2 



(4.30) 
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Again when we consider the function G-{t) we get a power law of t 



5- 



2T[iC2t 



2 e*" 



»2(l+tanh2(ei2)) 



(4.31) 



In this case the presence of the constriction affects the exponent of this 
correlation function and can change the behavior of the tunneling amplitude. 

The two limits of short and long times (with respect to to) are clearly 
visible in the numerical evaluation of the function S-{t,d) as shown in the 
Fig. 14.11 In the numerical calculation we have fixed the value of d and 612 and 
then varied the value of C2t. As it is seen from the Fig. 14. II the two limits we 
have discussed are reached for C2t ^ d and C2t <tC d. In the Fig. I4.1l we report 
the behavior of the functions ReS and ImS vs. the time for some values of 
d and fixed Ou- Notice that we have restricted our calculation to the case 
C2t/L2 <^ 1 and d/L <^ 1. The agreement of the calculated expressions and 
the approximated results ()4.27|1 and ()4.3()j) is very good. 

Hence from now on we approximate the S- {t, d) function as the sum of 
the long and short time behaviors 



S_{t, d) =Q{t - to)S^{t, d^O) 
+ 0(to-t)(5-(t,rf- 



00 



- A_ 



(4.32) 
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Figure 4.1: a) Plot of ReS-{t,d) vs. In{c2t/L2) for various values of d. 
Observe the two different regimes for C2t > d and C2t < d. The two slopes 
agree very well with the approximated result of Eq. ()4.24|) . We have chosen 
exp(6'i2) = 1.5275 in this calculation, b) Plot of ImS{t,d) vs. In(c2t/L2) 
for various values of d. We used the same parameters as a). The values for 
small and large 02^/^2 agree well with the expected results (see Eqs. fl4.27|l 
and ()4.H()|l ). The downward curvature at large times arises from the finite 
size of the system used in the numerical calculation and disappears in the 
limit of large system size. 



where the A_ function assures the continuity of S- at the point t = to. The 
two functions S'_(t, d ^ 0) and S_{t, d — > 00) are given by the limits for long 
and short times of the function ReS and ImS, respectively. 

With this approximation we have separated the calculation of the tun- 
neling correlation function G- in the calculation of the short and long times 
behaviors. As we have seen these behaviors are similar to the translationally 
invariant case. We then expect that the low energy behavior (which corre- 
sponds to the low bias voltage region) of the conductance will be dominated 
by the long time part of S-{t,d). Vice- versa, the response to a high bias 
voltage will be dominated by the short time behavior of S-{t,d). Within 
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this approximation the G-{t,d) reads 



G_{t) = Q{t - to)G'_,,^oW + e(to - t)G'_,,^oo(t) exp 



2e* 



-202 



(4.33) 



where the exponential factor in the second term in the right hand side stems 
from the presence of the constant A_. 

We have thus obtained the expression for the G- correlation function 
when the constriction is present. The next step is to evaluate the function 
gT{ujT)- In this case we must evaluate an integrals of the form 



dt e"^\5±ity 



(4.34) 



io 



In the Appendix ini we discuss some details about this calculation. Our result 
for the gxiuJT) function is 
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(4.35) 



where we have defined 



2e 



*2 



a = — - 
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(4.36) 



/5=-Ve-2^^(l + tanh2(0i2)) 



and r(zi, z-i) is the incomplete F function [2l|. Notice that in the limit Q\2 — > 
we have a = /? and we recover the result of the translationally invariant 
case. Notice also that this function is written as the sum of two distinct 
contributions which came from the long and short time regimes. However 
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the gT function does not show two distinct behaviors for small and large 
frequency as is shown in Fig. 14.21 where we plot the longitudinal resistance 
Rxx, defined in the Eq. ()3.42|1 . for various values of the mixing angle 6*2 inside 
the constriction. In fact the Fourier transformations of both the long and 
short time regimes have contributions in the whole range of the frequency 
values. 




(0^[a.u.] 

Figure 4.2: Plot of the resistance Rxx/\^\'^ given by Eq. ()3.42|1 with the 
gx calculated in Eq. ()4.35p for various values of 02 at fixed 9i = 0. The 
oscillations at large bias voltage becomes more and more pronounced with 
increasing 62- 



In Fig. 14.21 we plot Rxx{'^t) in the case that the inter-edge interaction 
is confined to the region of the constriction (i.e., we set 6'i = ^3 = and 
let 62 assume several different values). Experimentally, 62 can be increased 
by narrowing the constriction by the application of a gate potential. When 
^2 = there is no interaction and Rxx diverges as at low bias. 

This low-bias behavior does not change upon increasing 62 because the 
long time behavior is dominated by the exponent a which does not depend on 
^2- At larger bias voltage, however, the plot of Rxx shows oscillations, which 
become more pronounced with increasing 92. We can express the period of 
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these oscillations in terms of the physical parameters of the theory 



A ^ . h hci cosh 26i , , 

^Vt = ^ = -^—tt^ ■ 4.37 

e*to e*c(cosh26'2 

The frequency of the oscillations increases with increasing 62 as it is apparent 
in Fig. 

Making the conformal transformation ()4.20j) in Eqs. ()4.28|) and ()4.31|) . 
and substituting in Eqs. ()4.24p and ()4.33|) we obtain, after lengthy calcula- 
tions, the result for in the case of finite temperature 



\2tic2J V 2c2 J sinh°(7rTto) 

e*^T*o / a . ujt 1 

+ a-t^ (^a, 1; 1 + - - j^^^o 

pA'A + t: - 



(4.38) 



p-i^ V 2 27rT' 1 -e2-^*o 



where F is the hypergeometric function of four arguments (also indicated 



as 2-^1) and B the Euler beta function [21[. In the case 9i = 62 we have 
a = /3, the first and third term cancel against each other and we recover 
Wen's result. 

In Fig. I4.3r a-d) we plot the differential resistance Rxx vs. bias voltage 
for a system without inter-edge interaction (dashed line - ^2 = ^1 = 0) and 
with inter-edge interaction (solid line - 9i = 62 = 1) for different values 
of irdT/ci = 0.1,0.5,1, and 1.5. The non vanishing value of 62 within the 
constriction induces oscillations in the Rxx vs. ujt relation with the same 
period as in the zero temperature case. However, we now have a maximum 
at zero bias voltage and two minima at finite bias voltage. This behavior is 
due to the fact that the temperature introduces a new energy scale. When 
the e*V > ksT we are essentially in the zero temperature case and the 
resistance Rxx decreases with decreasing bias voltage (see Fig. 14. 2|) . But, 
when the e*V < ksT the resistance turns around and begins to increase, 
reaching a maximum at zero bias. This behavior implies the presence of two 
minima located at bias voltages of the order of magnitude of ksT/e*: these 
are clearly seen in Fig. 14.31 The finite value of Rxx at zero bias (independent 
of Vt to first order) indicates that the constriction is behaving like an ohmic 
resistor in this regime, even though the resistance is strongly temperature- 
dependent. 
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Figure 4.3: Plot of the differential resistance Rxx/\^\^ vs. ujt for a system 
with inter-edge interaction within the constriction (continuous line, 6'i = 0, 
O2 = 1) and without inter-edge interaction (dashed line, 61 = 62 = The 
four curves correspond to different temperatures: nTd/ci = 0.1 (a), 0.5 (b), 
1 (c), and 1.5 (d). 



The presence of a constriction adds another energy scale in the problem, 
associated with the inverse of the characteristic time to- For temperatures 
smaller than h/to the low bias behavior is dominated by the same expo- 
nent a (cf. Fig. 14.31 (a.b)) irrespective of whether the inter-edge interaction 
is present or not. When the temperature, instead, is greater than h/to the 
exponent f3, which depends on the strength of the interaction within the con- 
striction, controls the behavior of Rxx (cf. Fig. I4.3r c.d)). As a consequence 
the minima at finite bias are generally deeper and shift to lower voltages. 

The effect of the constriction depends quantitatively on both the inter- 
edge interaction parameter 6*2 and the temperature. To appreciate this we 
plot in Fig. I4.4l the differential resistance Rxx for different values of the inter- 
edge interaction and the temperature. More specifically we have plotted R^x 
without interactions {62 = 61 = 0) and with interactions within the constric- 
tion {01 = 0, 62 = 0.2) for ndT/ci = 0.5,1,5, and 10. We notice that the 
effect of the inter-edge interaction disappears at sufficiently low temperature, 
since it is always the long times exponent a that matters in that regime. The 
effect of the interaction shows up upon increasing the temperature above the 
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Figure 4.4: Plot of the differential resistance Rxx/\^\'^ vs. frequency with and 
without inter-edge interaction within the constriction. Solid line - = 0, 
62 = 0.2; Dashed line - 61 = 62 = 0. Temperatures are irdT/ci = 0.5 (a), 1 
(b), 5 (c), and 10 (d). 



crossover energy h/tQ-. the latter decreases with increasing 62. Such a trend 
is clearly seen by comparing Figs. l4.Hl and l01 
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Chapter 5 
Conclusions 



In this chapter we compare our results with the theory of Wen and the 
measurements made by a group in Pisa showing that our model can explain 
some of the experimental features that are missed by the previously developed 
theories. After that we would like to discuss a few future research lines we 
can follow to further investigate these systems. Finally, we will state our 
conclusions. 



5.1 Comparison with the experiments 

The group of F. Beltram, V. Pellegrini, and S. Roddaro in Pisa has performed 
a series of measurements of the tunneling current due to quasi-particles tun- 
neling between two edges belonging to the same Quantum Hall liquid [28^ . 
In such a system, it is believed that the tunneling particles are the bulk 
quasi-particles with fractional charge. The results of Wen and Kane 
and Fisher jl^, predict that the tunneling conductance G{T,V) must show 
a power law behavior both in voltage and in temperature, 

G(0,\/) ~ 1/"; G'(T,0)~T" (5.1) 

with the same exponent a determined by the filling factor 

a = 2u-2, (5.2) 

and, because in these theories u < 1, we have a < 0. The divergence 
of the conductance at very low bias is removed by the presence of a finite 
temperature, thus recovering a linear relation between the current and the 
voltage. Notice that the expected divergence is G(0, Vr) —oo for Vr —* 
then the effect of the temperature is to generate a zero bias maximum and 
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two minima at finite bias. In the Fig. 15.11 we report the experimental data 
for the measure of the tunnehng conductance [28|. The magnetic field is 
determined by fixing the filling factor to z/ = 1/3 while the temperature is 
varied from 900 mK to 30 mK. 



B=6T 




-3-2-10123 

V [mV] 

Figure 5.1: Plot of the differential conductance vs. the tunneling voltage 
for various temperatures at fixed magnetic field. The magnetic field and the 
experimental device fix = 1/3. Data from j2^. 



For high temperatures (900 — 400 mK), the central maximum is clearly 
seen. Also two deep minima are present. The general behavior is in quali- 
tative good agreement with the expected theoretical behaviors (see Figs. 14.31 
and 14. 4j) . In this range of temperatures is also evident an asymmetry for the 
change of the sign of the tunneling potential. The origin of this asymme- 
try, which becomes more and more evident lowering the temperature, is not 
completely understood. There is the possibility, as discussions with the ex- 
perimentalists have pointed out, that it can be an artifact of the experimental 
setup. 

For lower temperatures (400 — 30 mK) the structure of the response func- 
tion changes dramatically. The asymmetry becomes more evident, the central 
maximum disappears and a deep central minimum starts to develop. The 
origin of this central minimum is completely unknown and we expect that. 
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if it is a genuine experimental result, we must plug in some other physical 
insights in our model to explain it. 

For high temperatures and in the region of the central maximum, both 
our model and that of Wen can be compared with the experimental results. 
To do that we need to fit the experimental data with the theoretical curves 
to fix the free parameters: the translationally invariant interaction angle for 
Wen or the inside and outside mixing angles in our model. We choose to fit 
the data by using the least square minima technique. We have then looked 
at the minima of the sum of the squared difference of the experimental and 
theoretical curves to fix the parameters. 




-2 -1.5 -1 -0.5 0.5 1 1.5 



V [mV] 

Figure 5.2: Plot of the experimental and theoretical results. The filling 
fraction = 1/3 and the temperature is fixed at T = 500 mK. We have fixed 
the total amplitude of the theoretical curves because we do not known the 
value of the tunneling amplitude |r| and partially removed the asymmetry 
by moving the central maximum at Vr = 0. 

In Fig. 15.21 we report the curves for the longitudinal conductance at T = 
500 mK obtained after having determined the best fit parameters. From 
this plot is clearly seen that both the Wen's theory and our model recover 
the structure of the central maximum. However the Wen's theory does not 
recover neither the deep minima that are present at \V\ ~ 0.7 mV nor the 
rise of the conductance at larger voltage. 
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We must point out that both the models predict either an attractive 
interaction in the region outside the constriction or a larger than e* = e/3 
quasi-particle charge. In fact we use as fit parameters the exponents of the 
tunneling conductance that tie together the quasi-particle charge e* and the 
outside interaction parameters in products like (e*/e)^e~^^ and thus it is 
not possible to obtain from our fit direct measurements of the quasi-particle 
charge and of the mixing angle 9, separately. Similar fits with the other 
experimental curves in the range 900-400 mK have shown an increase of the 
term (e*/e)^e~^^ when lowering the temperature. This dependence is not 
understood and is not predicted neither by our model nor by the theory of 
Wen. On the other hand it is not possible to extract from the experimental 
data a measure of the oscillation period as determined by the Eq. ()4.37|) . 



5.2 Perspectives 

The results obtained so far allow us to consider some possible developments. 
In particular two arguments are attracting our attention. The first is the 
study of the tunneling in the edges of the Quantum Hall liquid in more 
detail. Indeed, when we have dealt with the problem of the tunneling into 
an edge, the problem of the fine structure of the edge was not considered 
at all. However has been pointed out by several authors that, even for the 
Laughlin state where v = 1/m with m an odd integer, the edge can undergo 
a reconstruction li^ 41, li^, li^ and acquire a complex structure. This 



reconstruction can affect the tunneling amplitudes by modifying the form of 
the quasi-particle creation and annihilation operators. One way to take into 
account this effect can be the inclusion in the tunneling Hamiltonian of other 
terms that describe the possibility that several species of quasi-particles can 
propagate into the edge. On this research line lies also the result of Eric 
Yang et al. The quasi-particle species that can undergo the tunneling 

between two edges depend on the separation width of the edges. In particular, 
if one considers the tunneling between two Laughlin states, when the edges 
are far-off the tunneling quasi-particles are the fractionally charged Laughlin 
quasi-particles. However, when the edges are close enough, other quasi- 
particles with different charges can do the tunneling. This research line can 
thus clarify the properties of the tunneling quasi-particles and their relation 
with the bulk quasi-particle. 

The other argument we want to explore is a more detailed study of the 
effect of the constriction on the transport properties of the Quantum Hall 
bar. In fact we have not yet considered the possibility that the constriction, 
behaving like a semi-opaque barrier, has some tunneling resonances that will 
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enhance the longitudinal conductance even when the magnetic field is not 
present. A possible way to attack this problem is to solve the quantum 
mechanical problem of the tunneling amplitudes of a constriction in a bidi- 
mensional system with zero magnetic field. We expect that these amplitudes 
depend on the energy of the impinging wave. When the magnetic field is not 
zero one can use the basis developed for the tunneling in the previous case to 
calculate the new transmission and refiection coefficients. This research can 
also open the way to a more accurate study of the dependence of the tun- 
neling amplitude F as a function of the physical parameters of the impinging 
quasi-particles. 

5.3 Conclusions 

The main subject of this thesis has been the study of the transport properties 
of a two dimensional electron gas in the presence of a high magnetic field. 
When the magnetic field is very strong the kinetic energy is quenched and 
the dynamics of the interacting electrons is produced by the commutation 
relations that appear among the density fiuctuations as a result of the pro- 
jection into the lowest Landau level. The phenomenology of the transport 
properties is shown in the Fig. 11.11 As we have discussed the transversal 
resistance R^y (also known as the Hall resistance) shows a series of plateaus 
at values an = I'e^ /h where v is either an integer or a rational number, while 
the longitudinal resistance Rxx shows an alternation of vanishing regions and 
5-like peaks. This effect takes the name of Quantum Hall effect due to its 
connection with the Classical Hall effect. 

To put our work in context, we have reviewed some of the theoretical 
approaches developed for the calculation of the response function of these 
systems. In this thesis we have chosen a hydrodynamical approach that 
allows us to treat z/ as a free parameter without restricting ourselves to the 
case of V being an integer or a rational number. This allows us to describe 
the dynamics of the quasi-particles as charge density fluctuations localized 
at the edges of the system. We have then developed the algebra to calculate 
the response functions by using these density fluctuations and its relations 
with the quasi-particle creation and annihilation operators. Our approach 
gives us the possibility to introduce, in a clear way, the interaction between 
the quasi-particles. We have shown in detail how our model can recover the 
linear relation between the fllling factor v and the transversal conductance. 

This research started as a collaboration with the experimental group of F. 
Beltram, V. Pellegrini, and S. Roddaro in Pisa. The focus of the experimental 
research was the development of devices to detect and use the quasi-particle 
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fractional charges. The experimental device used is reasonably well modelled 
by our idealized scheme (see Fig. 12. 3^ . In studying this problem we had 
to face with two distinct problems. The first has been the possibility of 
tunneling between the edges and the other has regarded the possibility that 
the presence of the constriction can affect the transport properties. 

We have thus decided to start by including the effect of the constriction 
directly in the Hamiltonian by considering the inter-edge interaction poten- 
tial in a piece-wise form. The model in such a case can be easily solved. We 
have shown that the presence of a finite size constriction does not affect the 
low energy transport properties thus recovering the ideal relation between 
the conductance and the filling fraction. This result can be physically un- 
derstood by observing that in the low energy limit, which corresponds to the 
long wavelength limit, the constriction becomes fully transparent and the 
transmission coefficient goes to 1. 

We have then discussed the possibility of tunneling between the edges 
inside the constriction. The presence of the constriction introduces a time 
scale, to = d/c. to is the time a charge density wave needs to travel through 
the constriction. This time introduces then a separation in the behaviors for 
long and short (compared with to) times. We have discussed the fact that 
for short times the effect of the constriction becomes negligible. The low 
energy behavior is mainly determined by the long times behavior and this in 
turn is modified by the presence of the constriction. We have compared our 
work with the previously developed theory of Wen, who considers the case of 
tunneling between the edges in a translationally invariant system 
Many differences have been pointed out. In the regime of low bias voltage 
both theories predict the presence of a large central maximum. Near this 
maximum two deep minima develop. The structure of these minima differs 
in the two theories: in particular we predict that these minima must be 
much deeper than in the Wen's theory. At larger bias voltages we predict 
the presence of small oscillations whose period 5V = h/e*to allows a direct 
measure of the traveling time and of the quasi-particle charge. 

We have then discussed the comparison of our result with some recent 
measurements of the tunneling conductance. We have shown that our model 
can explain some characteristics of the experimental results. Finally we have 
described possible future research lines. 
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Appendix A 



Proof of the commutation 
relations 



In this appendix we discuss the commutation relations for the 6pa operators 
and show how to obtain it directly from the Lowest Landau Level projection 
and the hydrodynamical approximation. We consider separately the two 
cases of sharp and smooth edge. 

A.l Sharp edge case 

Let us consider a single left edge. We imagine that the region y < is empty 
and there is a sharp wall at y = 0. Then we have the condition on the wave 
number A; < 0. Wc suppose that the length scale of variation of the confining 
potential in the y direction is smaller than /, hence we can neglect the states 
with k around 0"^ and approximate the number density with a 6 function 



Notice that the negative value of k are given by the condition y + kt^ = Q 
which identifies the maximum of the gaussian packet in the LLL and then 
the position of the particle. There is a subtle point in this definition. The 

function n{k) is discontinuous in A; = 0. We take n(0) = v/2. This as- 
sumption can be justified using the definition of the theta function such that 
6'(0) = 1/2. We define the transverse fiuctuation 



n{k)^u9{-k). 



(A.l) 



oo 




(A.2) 



— oo 
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so we have the commutation relation between this new operators 

oo 

[5p{x),5p{x')] = j dy dy' ^ (clcidh,m - c[c^dk,i^ 
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(A.3) 



- ^l{v)^,{v)^l{v')^n{v')\ . 
If we use the expression for the function (pk{^) in the LLL, in the given gauge, 
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we obtain 



oo 

/ 



dy (pl{x,y)iph{x,y) = e 



—i{k—h)x 



exp 



(A.5) 



and this allows us to calculate the commutator 

{k - hff 



ki^h 

X exp 
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exp 



2;2n 



qH 



1^0 lu d , 

The last expression has been obtained by observing that in the limit / — > 
the gaussian may be approximated by a delta function. Notice also that h 
is evaluated at — and q which corresponds to A; = 0. This leads to a factor 
1/2 according to the definition of n(0). 

In the presence of two edges the derivation is similar but in this case 
the values allowed for the wave vectors in the expressions for the density 
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fluctuations are different. Is such a case the Kronecker delta given by the 
commutator of the fermion operators cannot be fulfilled if the densities belong 
to different edges. Then the commutator is zero. Hence we obtain, if a,/3 G 
{L, R} are the indexes of the edge, 

[Sp^{x),5pfsix')] = -^^S{x - x'X^. (A.7) 

The appearance of the Pauli matrix is justified by the different forms of 
the na{k) functions. In fact for the number densities we have, if D is the 
distance between the edges, 

hence when we integrate over h they assume a different sign 

oo 

j ^ [uLih + q)- n,{h - q)] = (A.IO) 

— oo 
oo 

^[nn{h + q)-nn{h-q)] = ^. (A.ll) 

This completes the proof for the commutation relations in the case of a sharp 
edge. 

A. 2 Smooth edge case 

Next we will consider the case of the smooth edge. In this case we rewrite 
the equation ()2.19|1 in the form 

[5p(r), Sp{r')] = J2Hk) - n{h))4{r)Mr'Wi{r')Mr) (A.12) 

k,h 

and in the case of slowly varying density function we have 

n{k)c^n{h) + {k-h)dkn{k). (A. 13) 

The commutator of the density fluctuations reads 

[6p{r),6p{r')] = J^i^ - h)d,n{k)4{r)Mr')vK^')Vi{r) 

k,h 

=° -^dAn{k)6{y-y')6ix-x'). (A.14) 
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Now by remembering that n{k) = 27r/^po(2/) kl^ = y ^ dk = dy/P we 
have 

Mr), 5p{r')] = -ifd,{dyPo{yMx - x')5{y - y') (A.15) 
and integrating over the edge we finally obtain 

[5p^{x),5pp{x')] = f dy f dy'[5p{r),5p(v')] 

J a J 13 

= -il^dyal^pp5{x-x') (A. 16) 

where po(?/) is the equilibrium density and p is the bulk density. The appari- 
tion of the (T^ matrix is justified by the order of the limit of integration in 
the edge: We have 

^ dy dyp{y) = p{d) - p(0) = p, (A. 17) 

/ dy dyp{y) = p{R + d)- p{R) = -p (A.18) 

JR 

if d is the length over which the density goes from zero to the bulk value 
(supposed for simplicity equal for the two edges) . The final observation that 
f — 27r/^po completes the proof. 
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Appendix B 

Properties of the eigenvalues 
problem 



In this appendix we want to study some analytical properties of the equation 
()2.H4j) . First of all let us define the operators 

M„,^ = (B.l) 

oo 

Ho^^p = tJ- / dx d^Va,p{x,x')d:c'- (B.2) 

In J 

—oo 

With this definition we rewrite the equation of motion ()2.34p in the compact 
form 

uMi^ = Hip. (B.3) 

It is easy to see that H and M are hermitian operators and we request that 
H is positive definite (this will assure the stability of the physical system). 
Let us define the auxiliary function 

m = H^if (B.4) 

which is a solution of the equation 

= (h-^^MH-^) ^ = (B.5) 

if (y9 is a solution of (j2.34|l . Because M is an hermitian operator we have the 
results: 

1. the set of solutions forms a complete base of the Hilbert space, 
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2. the orthonormality condition is 

J2 / dx%,^{x)'i'm,a{^) = ^n,m, (B.6) 

3. the completeness relation 

n 

Because there is a one-to-one relation between ip and \i/ we have the 
following properties of the solutions of equation ()2.34j) : 

1. the solutions (f form a complete base of the Hilbert space, 

2. they are orthogonal with respect to the scalar product 



3. they satisfy the completeness relation 

-i'^UJnVnAA'friA^")^/^^^ = - x') • (B.9) 

n 

We obtain the relations reported in the text if we normalize the functions (pn 

as ifn/ A\^n\- 

Now we want discuss the degeneracy of the eigenvalues of the equation 

• If fnAA is a solution with given eigenvalue Un then the function 
VnaiA is also a solution with eigenvalue Lj_n = —oJn- 

• If ^mAA is a solution with given eigenvalue then the function 

pifrnAA is also a solution with eigenvalue —uJm- 

Then we have that if ^mAA is a solution then o'^^is'^mA''^^ ^^^^^ ^ solution 
with the same eigenvalue thus the solutions of problem ()2.34j) are doubly 
degenerate. 
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Appendix C 

Landauer-Biittiker transport 
formalism 



In this appendix we want to discuss briefly the Landauer-Biittiker formal- 
ism for the transport in one-dimensional system and its application to the 
Quantum Hall edge transport 4^ . 



We will discuss before the case of the Integer Quantum Hall phase where 
many edge states are present. Every such state is a possible channel for the 
current transport. In the clean case, when there is not tunneling between 
different edges, in a six probes geometry similar to the one showed in the 
Fig. 11.31 we can write the balance equations 



( h \ 
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-1 
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-1 



-1 \ 






1 / 



/ ^1 \ 



where we have assumed that the current in the i-th terminal is determined 
only by its potential and the potential of the contact directly "down-stream" 
with respect to the versus of propagation of the edge which leaves the termi- 
nal. Notice also that we have ordered differently the reservoir to have a direct 
comparison with the matrix conductance given by the Eq. ()3.41|) . We have 
also considered the possibility that + 1 propagation channels are present. 
To solve this set of equations one fixes the current of the terminal 2, 3, 5 
and 6 to zero then the total current is given by J = /i = —I^. We fix also the 
potential of the terminal 4 to zero. We obtain, by defining 0",^ = I /{Vi — Vj) 
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0^2,3 = Cr5,6 = 0, 

e2 (C.2) 

Cr2,5 = 0-2,6 = 0-3,5 = 0-3,6 = -J^{N + 1). 

It is possible to generalize such approach to the case of the FQHE by substi- 
tuting the number of possible channels + 1 with the filling fraction u. 

It is also possible to consider the effect of the tunneling between two 
edges. We model this tunneling by defining a transmission T and a reflection 
R probability. The reflection coefficient is given by the probability that a 
particle moving from 2 to 3 or from 5 to 6 tunnels to the other edge. We 
have 



/ h 






/ 


1 














-1 \ 




( ^1 \ 


h 











1 





-1 












h 








-1 





1 













V2 


h 














-T 


1 


-i? 







^3 


h 











-1 








1 









\h 


I 




V 








-R 





-T 


1 I 







(C.3) 



Notice that the gauge invariance and the current conservation require that 
R + T = 1. We can solve this system with the same requirements for the 
setup of the currents and voltages to obtain the conductance matrix in this 
case. In the limit R —>■ we recover the previous result. 

A more complex matrix can be obtained if one considers the reservoirs 
coupled as we have discussed in Eq. ()3.4H) . If we consider the non interacting 
case u = 1, V = from the Eq. ()3.4H1 we have 
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-1 + 5g43 
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(C.4) 



The gauge invariance and the conservation of the current fix the allowed 
values of the conductances Sg2i, Sg23, Sg^i, Sg^s: we have then 

R = Sg2i = SgA3 = -5g4i = -Sg23 = i lim ^^^^^ e~^^ (C.5) 

The presence of the constriction couples all the edges and we must solve the 
linear problem given by the Eq. ()3.41|) . In such case it is not possible to 
define a single reflection coefficient and express the solution only in terms of 
it. 
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Appendix D 
Evaluation of integrals 



In this appendix we provide a few details concerning the evaluation of the 
integral occurring in the calculation of the Fourier transform of the response 
function In the zero temperature case, we must evaluate an integral 

of the form 

oo 

(D.l) 

to 

where a is a positive real number. To do this we go in the complex plane 
of the variable t and consider, for positive frequency u, an integration path 
as the one shown in Fig. ID. II A specular path in the lower half plane must 



y 




Figure D.l: The integration path used to evaluate the integral ()D.1|) when 
the frequency is positive. A similar path, closed in the lower plane, is used 
in the case a; < 0. 

be used for negative frequency. We observe that the integrand function has 
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no poles in the complex half-plane of t with a positive real part. Hence the 
integral on the whole path is zero. The integral on the arc vanishes by letting 
the radius go to infinity. 
As a result we get 



dt{5 ± it)-"e^^* ^-^ -"-^ 



i(^l)-°cj"-^r(l - a, -iwto 



(D.2) 



to 



The case to = can be carried out by calculating the convolution product 
between the Fourier transform of the B(t) function and the integral 



oo 

I 



dt{5±it)-''e'^' =2sin(7ra)e*2"|,^| 



X r(l + a)(Tl) 



-l-a 



-1-a 



(D.3) 



obtained by cutting the complex plane along the imaginary axis, starting 
from t = ±i6 as is shown in the Fig. ID. 21 when the singularity is located in 
the upper plane and for positive frequency. Similar paths were used in the 
other cases. 




Figure D.2: The integration path used for the evaluation of the integral ()D.3|) 
when the frequency is positive. A specular path with respect to the real axis 
is used when the singularity is located in the lower plane. Similar paths were 
used in the other cases. 



In the finite temperature case, we need to calculate the integral 



dt e 



iujt 



to 



sinh(7rTt)'^ 



(D.4) 
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One can easily obtain the result reported in the text by means of the sub- 
stitution s = e"^'^^* which reduces the above integral to the definition of 
the hypergeometric F function of four arguments j2l|, . The integral with 
to = can be easily obtained in the limit to by using the corresponding 
limiting expression of the hypergeometric function F in terms of the Euler 
beta function. 
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